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ABSTRACT 


William  Walton  McDonald,  Master  of  Science  in  Engineering,  Arizona 
State  University,  May,  19 66.  Analysis  of  Stochastic  Networks ,  Major 
Professor:  A.  Alan  B.  Pritsker. 

The  purpose  of  this  research  was  to:  (l)  develop  a  general  purpose 
GERT  Simulation  Program,  (2)  investigate  resultant  distributions  of  GERT 
networks  containing  all  AND  nodes  to  determine  validity  of  the  PERT  nor¬ 
mality  assumptions,  and  (3)  investigate  analysis  of  AND  nodes  through 
analytical  and  simulation  rretnods  for  distribution  of  the  equivalent  time 
parameter. 

The  Simulation  Model  is  a  fast,  flexible,  user-orientud  computer 
program  requiring  only  one  input  card  for  each  activity  in  the  network. 
Five  probability  density  functions  are  available  for  use  in  describing 
the  distribution  of  activity  durations.  Program  output  includes  a  criti¬ 
cal-  .y  index  on  each  activity  and  node  analysis  on  specified  nodes ,  in¬ 
cluding  mean,  variance,  probability  of  realization  and  histograms  of  node 
realization  times.  Various  GERT  networks  were  simulated  to  provide 
examples  of  the  use  of  the  Simulation  Model  and  verify  results  through 
statistical  tests. 

Analysis  of  resultant  distributions  of  PERT  networks  verified 
the  PERT  normality  assumptions  utilizing  approximation  formulas  for  mean 
and  variance  of  activity  durations. 

The  merge  bias  correction  procedure  utilized  to  produce  better 
approximations  of  network  realization  times  was  found  not  applicable 
to  general  solution  of  GERT  networks  containing  AND  nodes. 
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CHAPTER  I 


BflRQDUenOH 

Stochastic  processes  are  defined  as  time  variant  and  probabilis¬ 
tic  processes.  Analysis  of  networks  with  stochastic  properties  has  been 
greatly  facilitated  through  the  application  of  GERT  (Graphical  Evaluation 
and  Review  Technique).^  This  procedure  has  provided  an  efficient  means 
of  analysis  and  communication  of  systems  problems.  A  complete  discus¬ 
sion  of  the  historical  development  and  theory  of  networks  may  be  found 
in  the  GERT  document. 


I.  STATEMENT  OF  THE  GOAL  AND  PR05LEH 

The  goal  is  to  develop  a  generalized  technique  for  analysis  of 
complex  systems  portrayed  by  stochastic  networks.  The  research  will 
concern  two  related  problem  areas. 

The  problem  is  to  develop  a  general  purpose  GERT  Simulator  to  allow 
analysis  of  GERT  networks  Through  the  technique  of  Simulation.  Analysis 
of  resultant  distribution  of  network  realization  will  be  facilitated  by 
this  simulation  program  to  assist  analytical  solution.  Specific  object¬ 
ives  will  be  to: 

1.  Accept  as  input  the  GERT  logic  rode  assignments  for  a  network 
while  allowing  a  choice  of  probability  density  functions  and  parameters 
for  activities  comprising  the  network. 

^A.  Alan  B.  Pritsker,  GERT :  Graphical  Evaluation  ar.d  Review 
Technique ,  The  RAJ.'D  Corporation,  RM-4973— lI  AS  A  (Santa  Monica,  April,  19 66) 
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2.  Provide  as  output,  a  criticality  statement  on  activities  lead¬ 
ing  to  realisation  of  the  network  and  analysis  of  several  designated  nodes 
through  histograms  of  node  realization  tines  and  statements  on  frequency 
of  node  realization. 

3.  Investigate  the  resultant  distribution  of  GERT  networks  with 
reference  to  theory  developed  in  the  literature,  specifically,  the  nor¬ 
mality  assumptions  pertaining  to  PERT  networks. 

inother  problem  is  to  investigate  analysis  of  AKD  logic  nodes 
through  the  use  of  analytical  and  simulation  techniques  for  probability 
of  node  realization  and  the  distribution  of  the  equivalent  time  parameter. 

II.  VALUE  GF  AH  AN SSEP. 

The  answer  is  of  value  in  that,  for  moderate  to  large  GERT  net¬ 
works,  the  computational  aspects  of  analytically  determining  the  distri¬ 
bution  of  network  realization  are  immense.  Simulation  cf  GERT  networks 
provide  a  fast,  efficient  means  for  obtaining  the  parameters  and  the  shape 
of  the  resultant  distribution.  The  criticality  of  an  event  contributing 
to  network  realization  provides  greater  insight  to  the  inne r-workings  of 
the  network.  The  complete  GERT  network  is  simulated  providing  probabil¬ 
istic  statements  about  the  terminating  nods(s) . 

The  basic  PERT  assumption  of  normality  of  project  completion  for 
large  networks  is  a  subject  of  controversy,  limiting  wide  acceptance  of 
this  technique.  Simulation  of  various  networks  under  fixed  and  probabil¬ 
istic  conditions  should  provide  an  answer  .ubstantuating  the  basic  assump¬ 
tion  or  rejecting  it,  and  pointing  out  the  need  for  further  research. 

Analysis  of  A!fl)  nodes  is  an  area  requiring  research.  The  value 
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associated  with  attainment  of  this  need  will  be  to  expand  the  scope  of 
'application  of  GER?  analysis  of  complex  systems  and  place  the  research 
endeavor  one  step  closer  to  the  goal  stated  above. 

m.  BACKGROUND  AND  LITERATURE  REVIEW 

In  the  past,  network  analysis  of  operational  systems  has  been 
concerned  with  project  scheduling  type  networks  (FSRT,  CFM,  etc.)  and 
signal  flow  graphs. 

The  PERT  type  analysis  treats  the  tine  of  each  activity  or  event 
as  a  constant  or  random  variable,  where  all  activities  contribute  to 
project  completion  or  realization.  The  technique  is  based  on  assumptions 
restricting  analysis  to  approximations  of  the  system  (project)  parameters 
and  confining  construction  of  networks  to  a  specific  form.  An  attempt  is 
made  to  deal  with  random  tine  variations  by  replacing  these  variables  with 
their  expected  times  and  reverting  to  a  deterministic  problem.  Based  on 
these  expected  times,  the  duration  of  the  longest  path  through  the  network 
is  utilized  in  making  estimates  of  project  completion  time.  The  distribu¬ 
tion  of  completion  time  ic  assumed  normal  by  invoking  the  central  limit 
theorem.^ 

Stating  that  the  above  estimate  of  project  completion  time  is 

3 

always  optimistic,  Fulkerson  proposed  an  improved  method  for  calculating 
an  approximation  to  the  expected  length  of  a  critical  path  in  a  PERT  net- 


opecial  Projects  Office,  Bureau  of  Naval  Weapons,  Department  of 
‘the  Navy,  PERT  Summary  Report .  Phase  1,  (Washington,  D.  C.,  July,  1958) 

■^D.  R,  Fulkerson,  Expected  Critical  Path  Lengths  in  PERT’  Networks , 
The  P.AND  Corporation,  RM-3C75-P8,  (Santa  Monica,  March,  19&2)  p,  1. 
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vcrk  by  considering  the  activity  duration  times  as  random  variables. 

4 

A  study  by  MacCrimnon  and  Ryavec  -  provided  an  evaluation  of  the 
errors  inherent  on  the  basic  assumptions  through  analysis  of  the  minimum 
end  Maximum  duration  time  of  parallel  network  paths.  Their  conclusion 
was  that  several  independently  parallel  paths  of  approximately  equal 
duration  was  found  to  give  th9  largest  error  in  PERT  calculated  mean  and 
variance .  Also,  cross  connected  paths  with  shared  activities  reduced  the 
errors  considerably.  The  presence  of  one  path  of  length  significantly 
longer  than  any  of  the  other  paths  reduced  the  errors  further. 

An  even  more  revealing  conclusion  was  drawn  by  Charnes,  Cooper  and 
Thompson'*  in  their  critical  path  analysis  of  PERT  networks,  utilizing  chane 
constrained  programming  methods.  Whenever  there  were  parallel  paths  that 
alternated  in  criticality  and  involved  sufficiently  different  times,  they 
found  that  multimodality  of  network  distribution  times  often  existed. 

Simulation  of  the  basic  PERT  network  was  first  performed  by  Van 
Slyke,^  through  application  of  the  Monte  Carlo  technique  to  obtain  more 
valid  statistics  of  system  behavior.  Activity  durations  were  treated  as 
random  variables  by  assigning  probability  density  functions  to  the  indi¬ 
vidual  activities.  A  criticality  index  of  the  relative  frequency  an 
activity  appeared  on  th3  critical  path  of  a  network  was  presented. 

U 

K.  R.  Mac  Crimmon  and  C.  A.  Ryavec,  An  Analytical  Study  of  the 
PERT  Assumptions ,  The  RAND  Corporation  (Santa  Monica,  December,  1962), 

5 

Charnes,  W.  W.  Cooper,  and  G.  L.  Thompson,  "Critical.  Path 
Analysis  Via  Chance  Constrained  and  Stochastic  Programming,"  0jo.  Res. , 

Vol.  12,  No.  3,  1964,  pp.  460-470, 

^R.  M.  Van  Slyks,  "Monte  Carlo  Methods  and  the  PERT  Problem,"  Op. 
Res.,  Vol,  11,  No.  5,  1963,  pp.  839-860. 
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Flowgraph  theory  has  been  applied  to  the  analysis  of  probabilistic 


systems,  utilizing  basic  properties  of  flowgraphs.  Realization  of  the 

flowgraph  network  results  from  consideration  of  all  the  transmittances  or 
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branches  of  the  graph. 

The  general  concepts  and  fundamentals  presented  in  GERT  provide  a 
convenient  basis  for  analysis  of  complex  systems  portrayed  by  stochastic 
networks.  The  research  endeavor  responsible  for  the  conceptualization  of 
CERT  is  a  continuing  process,  as  the  need  for  a  generalized  technique  to 
analyze  stochastic  networks  has  not  been  fully  realized. 

Conceptual  and  computational  problems  exist  in  GERT  analysis  of 
stochastic  networks.  Specifically,  no  general  method  of  analysis  exists 
for  the  AND  logic  node.  The  need  for  further  research  is  pointed  out  by 

g 

Pritsker,  while  presenting  concepts,  and  approaches  and  examples.  The 
use  of  the  AM)  node  in  GERT  is  the  same  as  in  PERT,  where  all  activities, 
or  branches  entering  the  node  must  be  realized  prior  to  the  realization  of 
a  branch  emanating  from  that  node. 

IV.  ORGANIZATION  OF  THE  REPORT 

The  following  chapters  of  this  report  are  organized  in  the  follow- 
itig  manner.  Chapter  II  presents  a  review  of  the  fundamentals  of  PERT  and 
GERT  to  establish  a  base  for  succeeding  chapters.  Chapter  III  describes 
the  GERT  Simulation  Model  with  discussion,  flowdiagrams  and  computer  list¬ 
ings.  Applications  of  the  Simulation  Model  follow  in  Chapter  IV,  covering 


n 

rPritsker,  op,  cit. ,  p.  16, 

g 

Ibid, ,  Appendix  B, 
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•problems  presented  in  the  GERT  Manual.  Also  included,  are  several 
variations  of  a  PERT  network  utilized  to  test  theory  developed  in  the 
literature.  Chapter  V  presents  analysis  of  AND  nodes  through  analytical 
and  simulation  treatments.  The  summary,  conclusions  and  recommendations, 
are  contained  in  Chapter  VI. 


CHAPTER  II 


’  REVIEW  OF  THEORY 

The  purpose  of  this  chapter  is  to  present  the  general  concepts  and 
theory  of  PERT  and  GERT  necessary  to  establish  a  foundation  for  the  mater¬ 
ial  in  the  following  chapters. 

I.  PROGRAM  EVALUATION  AND  REVIEW  TECHNIQUE 

The  Program  Evaluation  Review  Technique  (PERT)  is  a  method  of  plan¬ 
ning,  scheduling,  and  controlling  a  project  by  first  defining  all  signi¬ 
ficant  segments  and  activities  of  the  project  and  then  constructing  an 
interconnecting  network  of  nodes  and  arrows  depicting  the  various  time 
and  precedence  relationships  necessary  for  completion  of  the  project. 

The  time  duration  of  each  activity  is  estimated  by  three  parameters 
a,  m,  and  b,  where  a  is  the  earliest  (optimistic),  b  is  the  latest  (pessi¬ 
mistic),  ard  m  is  the  most  likely  completion  time. 

The  Beta  distribution  of  the  form, 

f  (t)  =K(t-af  (b.t  /  a  <  t  <  b 

*  0  otherwise 

(where  K,  ,  and  #  ,  are  functions  of  a,  m,  and  b)  is  assumed  to  repre¬ 
sent  the  distribution  of  activity  duration.  A  simple  linear  approximation 
for  activity  duration  was  derived  without  emperical  evidence  and  is  in  use 
today. ^ 

^Special  Projects  Office,  Bureau  of  Naval  Weapons,  Department  of 
the  Navy,  PERT  Summary  Report,  Phase  1,  (Washington,  D.  C. ,  July,  1953), 
Appendix  B. 
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=  expected  time  of  an  event 
=  (  ft+taotb  )  /  6 

VJIB(t  )  =  (  b-s  )2  /  36 

O 

Although  these  two  expressions  were  obtained  from  the  Beta  distri¬ 
bution,  the  usual  procedure  is  to  assume  normally  distributed  activity 
duration  times  with  the  parameters  given  above. 

The  time  at  which  a  node  is  realized  (achieved)  is  the  maximum  of 
the  durations  of  the  inwardly-directed  paths  to  that  node,  since  all  of 
the  activities  on  the  paths  directed  into  the  node  must  have  been  completed 
Outwardly-directed  activities  from  a  node  cannot  start  (be  released)  until 
the  node  is  realized.  The  project  duration  is  then  the  maximum  of  the 
elapsed  times  along  all  paths  from  the  origin  to  the  terminal  node.  The 
path(s)  determining  total  project  duration  is  designated  the  critical  path. 

Solution  of  PERT  networks  is  accomplished  by  deriving  the  critical 
paths  of  the  network  based  on  the  expected  times  of  activity  durations. 

The  stochastic  element  is  completely  ignored  through  solution  of  a  deter¬ 
ministic  network.  The  distribution  of  project  or  network  completion  times 
is  assumed  normal  using  the  central  limit  theorem.  Probability  statements 
concerning  project  completion  are  made  from  the  normal  distribution  v.ith 
mean  and  variance  derived  from  activities  on  the  critical  path, 

II.  GRAPHICAL  EVALUATION  AND  REVIEW  TECHNIQUE 

2 

The  Graphical  Evaluation  and  Review  Technique  (GERT)  is  a  tech¬ 
nique  for  the  analysis  of  generalized  networks  embodying  probabilistic  and 

2 

A,  Alan  B,  Pritsker,  GERT :  Graphical  Evaluation  and  Review 
Technique,  the  RAND  Corporation,  RM -497 3 -NAS A,  (Santa  Monica,  Aprii,  i960) 


$ 
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time  varying  characteristics.  The  general  term  associated  with  an 
activity  in  the  network  is  a  random  variable  and  a  probability  of  occur¬ 
rence  is  associated  with  the  realization  (acceptance)  of  the  activity. 
PERT  networks  are  in  a  sense  GERT  networks  where  the  probability  of  real¬ 
ization  of  an  activity  is  or.e.  Therefore,  all  activities  in  the  network 
Bust  be  transversed  for  network  completion  to  occur. 

GERT  networks  are  constructed  in  the  familiar  manner  as  an  inter¬ 


connecting  network  of  nodes  and  arrows  depicting  the  various  time  and 
precedence  relationships.  However,  the  form  of  the  nodes  indicating 
logical  relationships,  and  the  parameter  notation  for  the  activity  duration 
times  allow  GERT,  through  use  of  applicable  network  algebra,  to  achieve  a 
significant  development  in  generalizing  network  analysis. 

The  dual  parameter  of  time  and  probability  is  incorporated  into  one 
variable  allowing  use  of  applicable  network  algebra.  In  describing  the 
time  parameter,  GERT  utilizes  the  moment  generating  function  of  the  den¬ 
sity  function  (M.G.F.)  of  the  time  to  realize  an  activity.  The  product 
of  the  probability  of  realizing  a  given  activity  and  the  moment  generating 
function  is  named  ths  w- function,  and  is  defined  as  follows: 

wj(s)  =  pjVs1 


=  the  probability  of  realizing  activity  j 


the  moment  generating  function  of  the  time  to  realize 
the  branch  activity  j 


Flowgraph  theory  provides  a  convenient  basis  for  analysis  of 


stochastic  networks  us inf  the  w-function  from  GERT  and  the  logical  char¬ 
acteristics  of  flowgraphs.  Toe  logical  characteristics  utilized  are 
shown  in  Figure  1  for  the  series,  parallel,  and  feedback  cases  through 
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application  of  the  law  of  nodes  of  flowgraph  theory. 

Substitution  of  the  w-function  for  the  lower  case  a,  b,  and  c,  in 
Figure  1  provides  GERT  a  single  parameter  embodying  time  and  probability 
characteristics  while  permitting  utilization  of  the  logic  of  flowgraph 

thecry. 


The  equivalent  w-function  for  a  simple  network  as  discussed  provides 
information  such  as  li valent  probabilities  and  time,  through  the  follow¬ 
ing  relationships.  The  equivalent  probability  is: 

p^  =  w^(0)  »  by  setting  the  variable  s  =  0. 


The  M  G  F  of  the  equivalent  time  is  given  by: 

M^(s)  =  w^(s)  /  w^(0)  . 

The  relationships  of  network  type,  equivalent  w-function  and  M  G  F 

for  the  three  basic  networks  is  as  shown  in  Figure  2  , 

°  » 

Although  most  networks  can  be  shown  to  be  combinations  of  series 
and  parallel  equivalent  networks,  a  more  orderly  approach  to  network 
reduction  and  evaluation  is  through  use  of  the  topology  equation  of  flow- 
graph  theory.  The  topology  equation  for  closed  loops  is  as  follows: 

H(s)  =  1  -r  S(-1)™L  (m)  =  0  ,  m  =  1,  2,  3,  . 

c*l  * 

where  I^(m)  is  the  loop  product  of  m  non- touching  loops.  The  summation 
is  found  of  all  possible  combinations  of  these  loops.  Non-touching  loops 


are  defined  as  a  series  of  branches  or  activities  which  form  loops  and 
the  nodes  of  these  branches  are  non-touching  with  those  of  some  other 
loops.  Figure  3  illustrates  a  sample  network  employing  1st,  2nd,  and 


o.  28. 

* 
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-h— % 

a,  Series  B  =  aA  C  =  bB  =  abA 


c.  Feedback 


B  =  ap  +  bjB  =  k|/(  1  -  bp 
C  =  b2B  =  b^/Cl  -  bp 


FIGURE  1 


LOGICAL  CHA PA CT  E  R  T  £  T I CS  0?  FLO-VGRAPK  THEORY 


Network  Graphical  Paths  Equivalent  Equivalent 
Type  Representation  Function  ’V  M.  G.  F.  M  (s) 

V  v 


FIGURE  2 

NETWORK  REDUCTION  EMPLOYING  THE  TOPOLOGICAL  EQUATION 
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*  3rd  order  loops,  vhere  the  order  of  the  loop,  m,  is  as  above. 


( 

Ha  j- 
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V 

“4 

1 

i 

FIGURE  3 

SAMPLE  NETWORK  CONTAINING  FEEDBACK  LOOPS 

There  are  four  loops  in  Figure  3  of  order  1:  L^(l)  =  v^w ^ » 

1^(1)  =  L^(l)  =  w^wg;  I^(i)  =  w^wg;  and  L^(l)  =  w^w^.  An  example 

of  a  second  order  loop  is  L^(2)  =  L^(l)  L^Cl)  =  w^w£w^w5»  anci  a  third 
order  loop  would  consist  of  L^(3)  =  1^(1)  L^(l)  L^(l)  =  wpw2w5w6w9wi0* 

Th9  equivalent  w-function  for  a  network  or  section  of  the  network 
may  be  obtained  by  closing  the  network  such  as  is  shown  in  Figure  3  > 
with  the  dotted  activity  from  node  1  to  node  E  and  assigning  the  equiva¬ 
lent  w-function  for  the  network  to  this  activity,  1  /  w^(s).  A  character¬ 
istic  of  the  topological  equation  is  utilized,  where  for  any  closed  network, 

H(s)  =  0,  for  all  s. 

Alternatively,  w^(s)  could  have  been  obtained  directly  through 
application  of  Mason's  rule  or  the  loop  rule  for  open  paths. 

wt(s)  =  ^E?(path  i)  (l  +  !C(-l)m  (loops  of  order  m  not  touching  path  i)  ) 

(l+S  (-l)m  (loops  of  order  m)  ) 

m  • 

GERT  utilizes  node  shapes  to  indicate  logical  relationships  of 
activities  terminating  and  emanating  at  a  node.  Three  types  of  input  node 


13 

and  two  types  of  output  nodes  a.e  considered,  Figure  4  presents  «  summary 
of  the  GEKT  node  logic  for  input  and  output  types  cf  nodes.  Also  included, 
•re  the  six  possible  combinations  of  input-output  logic  nodes  utilized  in 
constructing  GER?  networks. 

The  discussion  in  this  chapter  was  intended  to  present  the  basic 
concepts  pertaining  to  PERT  and  GER7.  GEHT  is  a  general  network  evaluation 
technique  encompassing  the  principles  of  PERT  and  Flowgraphs.  ^hile  assur.p- 
ticns  are  inherent  in  the  solution  of  PERT  networks,  GERT  provides  a  flex¬ 
ible  analysis  framework  applicable  tc  a  vide  range  of  problems.  Additional 
treatment  of  the  theory  and  applications  may  be  foum  ^  n  the  GSRT  document. 


Bane 

Symbol 

Chirac is ris tic 

INFUT 

EXCLUSIVE-OR 

k3 

The  realization  of  any  branch  leading  into 
the  node  causes  the  node  to  be  realized; 
however,  ore  and  only  one  of  the  branches 
leading  into  this  node  can  be  realized  at 
a  given  time. 

INCLUSIVE-OR 

<J 

The  realization  of  any  brans  n  xeadir.g  into 
the  node  causes  the  node  to  be  realized. 

AND 

a 

The  node  will  be  realized  only  if  all  the 
branches  leading  into  the  node  are  realized 

OUTPUT 

DETERMINISTIC 

D 

All  branches  emanating  from  the  node  are 
taken  if  the  node  is  realized,  i.e.,  all 
branches  emanating  from  this  node  have  a  p- 
parameter  equal  to  1. 

PROBABILISTIC 

> 

Exactly  one  branch  emanating  from  zhe  node 
is  taken  if  the  node  is  realized. 

The  six  possible 

types  of  nodes 

ars  jo  JO  o  o  0  0 

FIGURE  4 


GERT  MODS  LOGIC 


CHAPTER  m 


THE  SIHOLATICH  dODEL 

The  general  i :  ncepts  of  GERT  provide  a  convenient  basis  for  analy¬ 
sis  of  stochastic  networks.  For  modest  to  large  size  networks  of  complex 
systems,  the  computational  aspects  involved  in  analytically  determining 
the  distribution  of  network  realization  tines  are  immense.  The  simulation 
model,  or  prog  raw,  was  designed  to  facilitate  analysis  of  stochastic  net¬ 
works  as  defined  in  G8RT,  and  allow  as  much  freedom  as  possil  le  in  depict¬ 
ing  networks  for  analysis. 

The  purpose  of  this  chapter  is  to  present  the  construction  and 
operation  of  the  simulation  model,  deferring  applications  to  a  later  chap¬ 
ter. 

The  simulation  model  can  be  separated  into  three  distinct  parts: 

(1)  Initialization,  where  the  variables  are  initialized  and  input  data  is 
processed,  (2)  network  simulation,  and  (3)  program  output.  Five  Fortran 
subprograms,  in  addition  to  the  main  programs  are  utilized  to  process 
input  data,  assist  in  the  network  simulation  and  output  resultant  statis¬ 
tics  of  the  simulation  trials. 

Output  from  the  program  includes:  (l)  A  criticality  index  on  each 

activity  as  the  relative  frequency  an  activity  appears  on  the  critical 

% 

path,  (2)  the  mean  and  variance  of  network  completion,  and  (3)  analysis 
of  specified  nodes,  consisting  of  mean,  variance,  histogram  of  node  real- 
izati  on  times  and  probability  of  realizing  the  node.  A  Chi  Square  Goodness 
of  Fit  Test  may  be  applied  to  the  distribution  of  node  realization  times 
and  printed  as  output. 
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this  chapter  is  divided  into  four  sections  to  facilitate  discussion 


of  the  aodel.  The  sections  are:  (l)  rules  for  establishing  the  network 
■odel,  (2)  subprograms,  (3)  method  of  simulation,  and  (4)  input-output. 
Method  of  presentation  is  discussion  followed  by  flovdiagram  and  computer 
listing  where  applicable. 

I.  RUIZS  FOR  ESTABLISHES  THE  NETWORK  MODEL 

The  following  rules  are  necessary  to  establish  restraints  on  the 
network  model  so  the  computer  program  will  function  as  intended. 

1.  There  must  be  one  originating  node  for  the  network,  number 
1,  although  one  or  more  terminating  nodes  are  allowed. 

2.  Mode  numbering  must  be  increasing  from  beginning  to  end  of  the 
network.  The  I  and  J  node  designations,  or  attributes,  of  an  activity 
will  be  numbered  such  that  J  >  I  in  the  general  case. 

3.  Mode  numbers  may  consist  of  any  three  digit  number,  with  the 
exception  of  099*  which  is  reserved  for  input  data  control. 

4.  Feedback  or  looping  of  an  activity  to  a  prior  point  in  the 
network  is  indicated  whenever  the  I  and  J  node  designations  of  an  activity 
are  such  that  J<  I. 

5.  One  exclusive  -  or  node  must  be  placed  in  the  feedback  loop 

when  a  one-activity  feedback,  or  self-loop,  condition  exists,  to  preclude 

% 

the  possibility  of  J  =  I  as  indicated  below. 


16 


6,  The  maximum  number  of  activities  and  nodes  is  confined  to  400 

each. 


II.  SUBPROGRAMS 

Pive  Fortran  subprograms  account  for  the  greater  part  of  data  man¬ 
ipulation  and  notably  simplify  the  main  program  logic.  They  are  Subroutine 
BETAD,  Subroutine  SET(JS),  Function  DIST,  Subroutine  HIST,  and  Subroutine 
CHI.  This  section  will  present  the  purpose,  theory  and  discussion  of  each 
of  these  subprograms.  Flowdiagrams  and  computer  listing  will  follow  each 
subprogram. 

Subroutine  BET AD 

The  purpose  of  Subroutine  BETAD  is  to  convert  the  PERT  a,  m,  and 
b  time  estimates  pertaining  to  each  activity  in  the  network  into  the  para¬ 
meters  c<  and  &  of  the  Beta  probability  distribution.  The  flow  diagram 
is  given  in  Figure  5  and  computer  listing  is  found  on  page  19. 

The  standardized  form  of  the  Beta  distribution  was  previously 
given  as, 

f  (x)  =  K  (x)  **  (1-x)  ^  0  <  x  <  1 

K,  c<  ,  $  =  constants  >  o 

=  0  elsewhere 

Equating  expected  value  and  variance  of  this  distribution  with  the 
standardized  approximations  based  on  the  a,  m  and  b  time  estimates  gives 
the  following 

E(x)  =  (*+  1)  /  (*  +  +  2) 

=  (  (a  +  4m  +  b)  /  6)  -  a 


b  -  a 


17 


V(x)  =  ( +  1)  (  a  +  l)/(oC  +  &  +  2)2  (<*  +  X  +  3) 

=  1/36. 

The  modal  value  m,  similarly  transformed  gives 
r  =  (m-a)  /  (b-a) 

=  C<  f  (o<  +  X  ) 

Elimination  of  from  the  equation  for  r  and  V(x)  results  in  a 
cubic  equation  as  given  below: 

<*3  +  (36r3  -  36r2  +  7r)o<  2  -  20r2c<  -  24r3  =  0 
Given  r  =  (m-a)  / (b-a),  the  cubic  in  c(  is  solved  through  the  Newton 
Rapheson  iterative  procedure^-  for  -  unique  value  of  e{)l, 

The  fact  that  there  is  only  one  positive  real  root  of  the  cubic 
equation  in  is  established  through  application  of  Descartes's  Rule  of 
Signs.2 

Successive  application  of  »  Newton  Rapheson  equation, 

^i+l  =  ^  -  f (  €^i)/f»( 

results  in  a  very  close  approximation  t^  the  actual  value  of  °<  sought. 

The  subroutine  is  designed  such  that  the  initial  value  of  oc',  or  c/0t  is 
UO  times  the  modal  value,  100  *  r.  A  maximum  of  100  iterations  is  allowed 
to  achieve  an  accuracy  of  at  least  two  decimal  places,  f(  °<  ^) ^  1/100. 

Comparison  of  calculated  values  of  c<  and  5  with  tabulated  values 

^Kaiser  S.  lunz,  Numerical  Analysis  (New  York:  McGraw-Hill  Book 
Company,  Ir.c,,  1957 )>  p.18. 

2Ibld. ,  p.  11, 


FIG  cT  5 

FLO'V  BIAGFAF  OF  SUBPOUT  HI  bztad 
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FTN  1,4  0* TF  03/30/^?  AT  105^25  ?AGP  NO.  1 

SJ^OUTI-VE  irin(^iLPiA,SH?U) 

D!*«£mS!3V  ACID) 

sta-u  iThWAiiov  s -zyr.'  coffhcifnt 

NCOJvT  :<! 

ALPHt  =  19-?. 

Ns  3 
M=1 

5  C  0  ^  T  i  ;  j  E 
GO  TO  Cc,7.4,9)»l 

6  A(i)  =  3f;. 

A<?>=-A<1> 

A<3)*>. 

A  {  4  )  s  C  . 

B=R 
X  =  A  C 1  > 

Ms? 

GO  TO  10 

7  STOREY 

12  A ( t ) s 1 . 

A<2)sST3Rfc 
A<3! 

A(4)s-?4,*?**3 
X=  A ( 1  ) 

BsALPHA 
Ns  3 
Ms  3 

GO  TO  10 

8  Fsx 

IF(F-.P1)11»11.13 

13  A(l)s3.*A(l) 

A  ( ? )  s  2 .  •  A  ( 2 ) 

A<3)s-4<  5) 

Xs A ( i  ) 

Ns? 

Ms  4 

10  00  2n  ^ s  1.  #  %> 

20  X sy *'>♦:,(  v  +  i  ) 

GO  TO  5 

9  C  0  M  T  l  M  *  i  £ 

FPsX 

ALPHAS  ALPh,A-P/FP 
IF  (MTO*/ viT-1  J  1)14.14*50 

14  NC^u'.Ts  VCDJ‘!T->1 

GO  TO  1 2 

11  GAM'lisAUPH*»(l.-  )/B 
RETURN 

50  P^I'aT  u  1 » F ,  A l  ^ 

51  F0^‘*.,'T  c  J  i  IS.  =  l2.A»in  H>-  ITh  ai.sp  A  *  r  1  0  .  4  ) 

CALL  hIT  . 

F\n 
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given  by  Mac  Crimmon  and  Ryavec,  where  V(x)  =  1/36,  indicates  very 
close  agreement,  as  shown  below. 

BETAD  REFERENCE 


Mode 

Variance 

oC 

Jf 

* 

X 

1/4 

1/36 

1.22 

3.66 

1.21 

3.63 

1/3 

i/36 

1.82 

3.64 

1.82 

3.64 

1/2 

1/36 

3.00 

3.00 

3.00 

3.00 

Subroutine  SBT(JS) 

k 

Subroutine  SET(JS)  serves  the  purpose  as  master  time  file  for 
ordering  activities  filed  in  the  array  SETS  (M,  K)  as  they  are  released 
and  scheduled  for  completion  in  the  network  simulation  cycle.  Through 
the  ordering  process,  Subroutine  SET(JS)  maintains  the  system  clock 
pertaining  to  simulation  of  the  network.  The  flow  diagram  is  given  in 
Figure  6  and  the  computer  listing  on  page  22  through  24 

As  each  activity  is  released  and  scheduled  in  the  simulation  cycle, 
it  is  filed  in  the  array  SETS  (M,  K),  and  Subroutine  SET(l)  is  called. 

The  activities  in  SETS  (M,  K)  are  then  ordered  by  a  successor  -  prede¬ 
cessor  arrangement  based  on  column  KRANK  (l),  the  column  containing  the 
scheduled  time  of  each  activity.  The  three  attributes  of  an  acticity 
stored  in  SETS  (M,  K)  are  the  I  and  J  node  designations,  and  tQ,  the 
scheduled  tirr  of  completion. 

This  subroutine  also  provides  a  marker,  MIC  (l),  signifying  the 

^K.  R,  MacCrimmon  and  C,  A.  Ryavec,  An  Analytical  Study  of  the  PERT 
Assumptions,  The  RAND  Corpciation,  RM-3408-PR,  (Santa  Monica,  December, 
1962),  p.  44. 

4 

Don  Deutsch  and  Philip  Wolfe,  '*A  Revision  to  GASP,  A  Simulation 
Programming  Language”,  (Unpublished  tern  paper,  IE  477g,  Arizona  State 
University,  May,  19^5). 
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FIGUPE  6 
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must  INITIALIZE  KNAT<CJ5) 

IF(INITS-1)27, 23,27 
2*  EOl’s  7/77. 

I0r20 

NOO=l 

H=3 

EOF  =  3338. 

WLF  =  0999. 

MX  =  I’Hi 
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RETURN 
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DO  25  K  =  l,IO 
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row  containing  the  activity  bearing  the  lowest  scheduled  time  and  next 
for  exit»  based  on  the  scheduling  rule  FIFO.  The  marker  KFA,  or  next 
row  for  assignment,  is  also  provided.  Both  marker  variables  are  deter¬ 
mined  within  Subroutines  SET(l)  through  the  successor-predecessor  arrange¬ 
ment  mentioned  above. 

A  more  complete  description  on  the  operation  and  use  of  Subroutine 
SET(JS)  may  be  found  in  the  referenced  paper  by  Deutsch  and  Wolfe. 

Function  DIST 

The  function  subprogram  DIST  is  a  prime  part  of  the  computer  pro¬ 
gram  in  that  the  variable  DIST  takes  on  a  time  value  as  a  random  variable 
of  the  required  probability  function  of  an  activity.  The  parameters  of 
the  population  distribution  are  defined  *d.th  each  call  of  the  subprogram. 
The  flow  diagram  of  this  subprogram  is  given  as  Figure  9  on  page  29, 
and  the  computer  listing  is  on  page  32. 

A  brief  review  of  sampling  theory  and  use  of  Monte  Carlo  methods 
is  offered  to  facilitate  description  of  the  probability  distributions 
available  for  use  in  the  computer  program. 

Simulated  sampling  involves  replacing  the  actual  universe  of 
items  by  an  assumed  theoretical  probability  distribution,  and  then  samp¬ 
ling  from  this  theoretical  population  through  use  of  random  numbers. 

To  draw  an  item  at  random  from  a  universe  described  by  the  proba¬ 
bility  density  f (x) ,  one  derives  and  plots  the  cumulative  density 
function  F(x) , 

Y  =  F(x)  =  /  f(u)  d  u 

as  described  in  Figure  7« 
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FIGURE  ? 

EXAMPLE  OF  DRAWING  AH  INDEPENDENT  SAMPLE 
FROM  A  CUMULATE  DISTRIBUTION  FUNCTION,  F(x) 

A  random  number,  having  equal  likelihood  of  lying  anywhere  in  the 
interval  (0,  l)  is  selected  as  the  entry  on  the  y  axis  and  projected  to 
an  intersection  with  the  curve  y  =  F(x).  The  value  x  corresponding  to 
th6  point  of  intersection  is  determined  to  be  the  independent  sample 
value  of  x  from  the  distribution. 

With  this  background  in  mind,  a  brief  description  of  the  proba¬ 
bility  distributions  available  follows. 

1.  Negative  Exponential.  Distribution 


y  -  f(x)  =  (l/u)e  ”x^u  '  x,  u  >  0 

=  0  elsewhere 

F()()  =  ex^  x,  u  ;>  0 

=  0  elsewhere 


x  =  DIST  =  -u  *  In  (RAND(RN)  ), 

where  RAiJD(RN)  is  a  uniformly  distributed  random  variable 
and  satisfies  the  relationships:  0  ^  RAND(RN)  1,  and 
u  is  the  mean  of  the  exponential  distribution. 

2.  Normal  Distribution 

f(z)  -  (l/  f~2T f  )  exp  (-z^/'2), 
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where  %  -  (x  -  a)  /  Cf  .  Using  an  approximation  to  the 
normal',  random  normal  deviates  can  be  obtained  from, 

%  —  (  -2  ln(RAKD(R!0  )  )*  cos  (  2  RAKD(RN)  ). 

Thus  z  is  a  normally  distributed  random  variable  with  a  mean 
of  0  and  standard  deviation  equal  to  1.  Using  the  trans¬ 
formation, 

DIST  =  x  =  cf  z  -  u,  a  normally  distributed  variable  with 
a  mean  of  u  and  standard  deviation  of  <f  is  obtained.  The 
ratio  of  3  /  u  should  be  less  than  1/3  to  make  the  possibil¬ 
ity  of  negative  numbers  small. 

3.  Discrete  Rectangular  Distribution 

The  discrete  rectangular  distribution  takes  on  integer  values 
between  A  and  B  with  equal  probability,  thus; 

DIST  =  INTF  (A+(B  -  A)*RAND(RN)  ),  where  A  and  B  are  derived 
from  the  relationships  u  =  (A  +  B)  /  2  and  =  (B  -  A)^/l2. 

4.  Constant 

The  time  duration  of  the  activity  is  defined  as  constant  and 
equal  to  the  value  x. 

5.  Beta  Distribution 

As  described  previously,  the  Beta  distribution  of  the  form, 
y  =  f  (x)  =  K  x*  (l-x)r 

a  <  x<  b 

=  0  elsewhere 

^Claude  McMillan,  and  Richard  F.  Gonzalez,  Systems  Analysis . 
(Homewood,  Illinois:  Richard  D.  Irwin,  Inc.)  p.  157. 
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has,  as  its  cumulative  distribution,  an  indefinite  in¬ 
tegral.  Tabulated  values  of  the  Incomplete  Beta  Function 
for  limited  values  of  and  X  may  be  found  in  Harter.^ 

The  method  utilized  to  draw  independent  samples  from  the  Beta 

*7 

probability  distribution  is  the  rejection  technique  described  by  Kahn' , 
which  differs  from  the  sampling  procedures  utilized  for  the  other  proba¬ 
bility  distributions. 

The  basic  procedure  is  as  follows  with  reference  to  Figure  8  .  A 
point  is  chosen  uniformly  at  random  from  the  rectangle  with  base  of  length 
b  and  height  M.  If  the  point  falls  below  the  graph  of  f(x),  accept  the 
abscissa  as  in  independent  sample  value  of  the  distribution.  If  not, 
reject  it  and  try  again. 

Restating  the  procedure  in  step  form,  and  again  referring  to  the 
graph  in  Figure  8. 

1.  Obtain  two  random  numbers,  and  R£. 

2.  If  f(a  +  b  R2)  /  M,  let  x  =  a  +  b  R2. 

3.  If  R1>  f  (a  +  b  Rg)  /  M,  pick  two  new  random  numbers,  R^  and 
Rg,  and  try  again, 

A  programming  aspect  of  this  sampling  procedure  is  noted  in  the 
seemingly  low  efficiency  of  selecting  samples  from  *’e  distribution. 

Should  the  range  b  be  very  wide,  and  the  value  of  M  very  small,  the 


H.  Ieon  Harter,  New  Tables  of  the  Incomplete  Gamma-Function  Ratio 
and  of  Percentage  Points  of  the  Chi-Square  and  Beta  Distribution.  Aerospace 
Research  Laboratories,  Office  of  Aerospace  Research,  United  States  Air 
Force  (Washington:  Government  Printing  Office,  1964),  table  3. 

7 

Herman  Kahn,  Applications  of  Monte  Carlo .  The  RAND  Corporation, 
A'£CU-3259,  (Santa  Monica,  April,  195-0,  p.  10. 
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probability  of  acceptance  of  the  variable  as  an  independent  sanple 
from  the  probability  distribution  would  indeed  be  snail,  since  is  a 
rando~  number,  0  4.  R^  ^  1.  Many  repetitions  of  this  rejection  procedure 
for  eac..  acceptable  sanple  could  occur. 

A  significant  improvement  in  the  efficiency  of  sampling  through 
use  of  the  rejection  technique  is  realized  by  restricting  the  random 
variable  to  the  range  0  <  R^<C  f(M),  where  M  is  the  modal  value  of  the 
distribution.  The  rectangle  Kxb  shown  in  Figure  8  is  then  the  actual 
space  we  utilized. 

Although  the  rejection  technique  could  be  utilized  in  sampling 
from  the  other  distributions  described  in  this  section,  the  more  direct 
approximations  utilized  require  much  less  time  in  computation  and  are 
therefore  utilized. 

A  sample  network  consisting  of  one  activity  defined  by  the  PERT 
time  estimates,  a,  m,  and  b,  of  the  Beta  probability  distribution  was 
simulated  to  test  adequacy  of  the  sampling  procedure  outlined  above.  The 
GERT  Simulation  Program  of  this  report  was  utilized  in  simulating  the  net¬ 
work. 

The  test  was  to  determine  if  symmetrical  time  estimates  would  pro¬ 
duce  a  reasonable  approximation  to  the  normal  distribution  as.  determined 

0 

by  the  Chi-Square  Goodness  of  Fit  Test. 

The  hypothesis  that  the  sample  distribution  is  normally  distributed 
about  its  mean  was  not  rejected  at  the  5/°  confidence  level  since  the  calcu¬ 
lated  X^  =  19.3^  is  much  less  than  the  critical  level  of  ( . 95 a 17 )  -  27.6. 
Subroutine  HTST 


Subroutine  HIST  accomolishes  the  task  of  converting  the  01 


:out 


I 


Peturn 


FIGUP2  9 

FLO?/  DIAGRAM  OF  FUNCTION  DIST 


FIGUP3  9 

FLOW  DIAGRAM  0?  FUNCTION  DIST 
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FTN  1.4  DATE  P3/3T/66  AT  105*23  PAGE  NO 

FJMCT l  ON  Dl$T(x»xl ,x?,X5» *4) 

COMMON  RM 

COMMON  STATSaCOP.O 

COMilON  SETS  ( « 0  •  ? )  «  l  1  T5 « 9JT  ,  H R£  MK  (  j  ) ,  |  NN  <  1 )  •  *LC  ( 1 )  •  Mr  A 

Nsy2 

IF  {4)1P#1C»3.? 

10  CO-TI  -Mt 
DlST=0.0 
RETURN 

12  GO  TO  (1»2»3*4»5»6*7*B»9»15  )#N 

1  nisT=-y*ionp  (RAMLHR  j) ) 

RETURN 

2  RlrRA\n<P\') 

R2  =  RAN-n(o^) 

V= < -2. **LOOF < Rl )> **0, 5 *COSF( 6.283*92) 

D|ST=V*SCR1F(X1)*X 

RETURN 

3  |sXl*(y-Xl) *R4NtHRv i 
D I S  T  =  I 

.RETURN 

5  CONTINUE 

IF(X)17,?1,17 

17  I F ( X 1 ) 1 8 » 1 y , 1 8 

18  CONTINUE 
XMOOEsV/tX+Xl) 

XMOOFsV  10p=**X*<l.-<MQi)c)**Xl 

19  Tt  R AM0 ( R*  ) 

Y  s  °  A  v  0  (  R '  ) 

Y  s  Y  *  X  M  0  U  P 

8E  T  A 5  T  *  *  X  *  ( 1 ,  •  T  )  *  *  X 1. 

IF{JJLTA-Y)19,2n,?(» 

20  niSr*X^*T*(X4-X.$) 

RETURN 

'21  D I  ST 
P  c  T  U  3  N 

6  no  TO  A 

7  GO  TO  4 

8  GO  TO  4 

9  GO  in  d 

11  GO  TO  4 

4  D I S  7  =  X 

RETURN 

14  FORMAT  (6"7:.<0  '!//> 

ENn 


statistics  into  a  standard  histogram  configuration  of  twenty  cells.  This 
subroutine  requires  only  identification  of  the  storage  array  and  number 
of  records  contained  in  the  array  as  input.  This  is  a  general  purpose 
subroutine  with  internal  decision  rules  for  setting,- up  cell  size  and 
location  for  the  most  convenient  grouping  of  data  and  for  output.  Figure 
10  contains  the  flow  diagram  of  Subroutines  HIST  and  the  computer  listing 
is  located  on  page  35* 

The  data  array  is  reviewed  to  obtain  expected  time,  variance ,  high 
time,  low  time,  and  range.  Based  on  these  values,  the  histogram  is  set 
up  such  that,  if  the  range  is  less  than  20,  the  first  cell  will  be  the 
lowest  data  value  truncated  to  an  integer  and  cell  increments  will  be  of 
size  one.  Should  the  range  be  greater  than  20,  the  midpoint  of  the  histo¬ 
gram  is  the  expected  value,  and  cell  increments  will  be  in  integer  steps 
of  RANGE/20  +  J. 

Output  statistics  are  expected  time,  variance,  high  time,  low  time, 
range  and  histogram. 

Subroutine  CHI 

¥ 

Subroutine  CHI  is  an  auxiliary  routine  used  to  measure  the  dis¬ 
crepancies  or  fit  between  the  simulation  output  distributions  and  a  theo¬ 
retical  probability  distribution.  Provision  is  made  only  for  .comparison 
of  symmetrical  distributions.  Input  data  is  required  as  described  in  a 
later  section  of  this  chapter. 

This  subroutine  is  utilized  only  at  the  analysts  discretion  by 
placing  the  subroutine  call  card  as  shewn, in  the  main  program  listing  on 
page  45, in  the  histogram  output  section  following  the  CALL  HIST  statement. 
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MAIN  PROGRAM 


FIGURE  TO 

FLO*.'  riAGP  .. :  c?  su3"cu:'i;:~  hist 
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FTN  1,4  DATE  03/30/66  AT  105823  PAGE  *3.  1 

SUBROUTINE  Hi  ST  (NT,  NS  pi) 

COMMON  RK 

COMMON  STATS< ipco,  < ) 

C0MV.ON  SfcTS(*0»5>#  !•*  I  T$, OUT , KRA  MK  ( 1  ) ,  JMN<  1  >  # »*LC<  1  > ,  MF A 
DlME^'SIDv  iIST3(?0) 

DO  14  !sl,?f 
14  HlSTGCDsn.u 
BlG=-9. 0E5P 
SMALl.s9.nP5u 
EXPTsO.O 
VXPT=0.0 
X NS  I Ms MS IK 
DO  20  fst.^SIi 
XsSTATS(I,^T) 

EXPTsEX?T+X 

vxpTsvy-n+x**/ 

BIG*MAxlFmtS,X) 

20  SNALL=‘M‘ IF (SMALL. X) 

EXP Tstv;>TA '-SIX 
VXPT  =  VXPT/x\SI'*-GXPr*EXpT 
PRINT  1 1 7 , E  x P T 
PRINT  llP.VXPf 
C 

RANGE*0 IG-SMAuL 
lF(RAf!0=)5n,na,fli 

80  PRINT  113 
RETURN 

81  lF(RMOc-2r.)86,54.?7 
86  NSTZFsi 

SI7E=1.0 

NSmsSnali  : 

NBnsNS‘14.19 

ASMAl.L*NSM 
GO  TO  *3 

37  S|7F  =  R/\f-F/?3, 

NSUE*M7E 
H3IZF=517E*.5 
!F(.MSl?  =  -N$I7E)d9, 87.80 
8R  N3lZFs'<St2F 
89  KlnPrR>9T 

NSMs’l  I  n3-9*  iS  I  Zc 
IF ( NSM ) 90  0 » °n i , 901 

900  NS**s n 

901  ASULLsNSM 

M3G  =  VS  >19* '.SIZE 
S I 7~s  NS  I 7R 

83  OD  8?  '=i.V5H 
X  s  S  T  A  T  S  (  \ ,  V  r  > 

IF(X-AN4ALL )85. 34, A  ; 

84  Ls(X-A?MALLJ/SIZE*1  . 

|F{L-2*>»2,4?.«3? 

8  9  ?  L  =  5  ’j 

GO  m  -2 

85  t  =  1 

8?  H  I  ^  T G  ( iw  )  s  H I  -1 T G  ( L  )  ♦  1  .0 

PR!  H  U^, (  I  ,  I=-MSvf‘.aGt  i?i?p) 


# 
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PRINT  11Q,<nISTo(L>,Ls1#?3> 

PRINT  Ilfc.SIR.SMALURAVGr 
RETURN 

113  r0SiAT(//2<>X.2?N^!ST0G^Ax?  IS  CANCELLED) 

114  FO«UT«rU) 

115  FORMAT  t  /?0 1 4  ) 

116  FORIATr5x#3Hpl3»7X,5HS^ALL»5X.5HRA.vGE/3Fl0.3) 

117  F0R:|AT(//2rx,4CH=X3rCTrD  THE  OF  COMPLETION  THIS  NODE  IS, Flo 
lift  F0«MAT</. 2r*, ^VARIANCE  TIME  OF  COMPLETION  THIS  NODE  iS.FlC 
119  FORMAT ( /20r4 . 3 ) 

.END 
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A  measure  of  the  discrepancies  existing  between  observed  and 

expected  frequencies  is  given  by  the  X^  statistic  where  the  data  is  grouped 

into  cells  or  events  of  the  total  sample  population.  For  Large  N,  or 

2 

sample  size,  the  sampling  distribution  X  closely  approximates  the  Chi 

2 

Square  probability  distribution.  .  X  is  given  as, 

X2  =  £  (0.  -  e,)^  /  e.,  with  k-3  degrees  of  freedom. 

jit  J  J  J 

Each  summation  represents  the  observed,  0^,  and  expected,  e^, 

frequencies  of  occurrence  within  a  cell,  where  there  are  k  cells. 

Expanding  and  simplifying  the  above  expression  where  N  =£  0 ,  =  £  e . 

j*f  J  J 

gives  the  following  expression J 

x2  =  (  So.2/..  )  -N. 

J*l  J  J 

This  subroutine  requires  the  location  of  the  data  storage  array, 
the  number  of  records  contained  in  the  array,  and  whether  or  not  to  read 
input-data  describing  the  theoretical  probability  density  function. 

The  data  in  the  storage  array  is  reviewed  for  calculation  of  the 
expected  time  and  variance  and  then  transformed  through  the  normalizing 
equation  z  =  (x  -  E(x)  )  /  c f  *  .  As  the  data  is  normalized,  it  is 

placed  in  a  20  cell  histogram  with  center,  z  =  0;  variance,  V  (z)  =  1; 
and  increments  of  3/10.  Therefore,  the  end  cells  of  the  histogram  repre¬ 
sent  +_  3d  •  The  grouped  histogram  cells  are  utilized  in  calculating 
2 

the  value  of  X  . 

2 

Output  of  the  subroutine  includes  the  value  of  X  calculated, 
expected  time,  variance,  and  the  normalized  histogram  of  data. 

Tests  of  Hypothesis  may  be  made  at  the  significance  level  desired, 


with  1?  degrees  of  freedom,  based  on  the  tabulated  values  of  the  Chi 
Square  probability  density  function. 


MAIN  PROGRA! 


SUBROUTINE  CHI 
Given  storage  loca¬ 
tion  and  size 


Review  data  for  u,<T 


Normalize  data,  place 
in  histogram-3<Z<+2 
with  center  at  zerc 


Read  input  data 
defining  theoret¬ 
ical  prob  function 


Compare  expected 
With  observed 
X2=2(o//eJ  - 


Print  X  .  mean.  Va 


Return 


FIGURE  11 

FLOW  DIAGRAM  OF  SUBROUTINE  CHI 
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SUBROUTINE  CHI <NT,NSM,NN» 

COMMON  RN 

COMMON  STATSUOOOO.l) 

DIMENSION  HIST (20 )  #  l)ATA(  20 ) 

DO  2  1=1*20 
2  HIST ( I  )  =  0 .0 
VARsQ.O 
SUM=0.0 
XCHI=0.0 
DO  10  1=1, NSM 
XssTATSU*NT) 

SUMsSUM+X 

10  VA.DsVAR*X*X 
XNSM=NSM 
ET IMF=SUM/XNSM 

var=var/xnsm*etime*etime 

SIG=SQRTF(VAR) 

DO  20  1=1. NSM 
XsSTaTSU.NT) 

ICOOE=1.0 
Z*'X»ETIKE>/SIG 
IE  <Z)5»7f 7 
5  ICODE=0 , 0 

7  NsaBSF (Z)*10./3«*l. 

IF(N-10>12,12,11 

11  N=10 

12  L=N*10* ICODE 
IF (L«10)8,8*9 

8  L=11»l 

9  H!STILl*HlSTCl>M. 

20  CONTINUE 

IF(NN-l)22, 22*41 
22  DO  40  1=1,10 
READ  66.DATAU) 

•  DATA(21-I>  =  DATAU) 

40  CONTINUE 

41  XCH!=0.0 
M=0 

DO  60  1=1*20 

XCHI=(HIST( I>**2)/tOATA(M*I)*XNSH)+XCHI 
60  CONTINUE 

XCwJsXCHI-XNSM 
PRINT  67.XCHI 
PRINT  70 , ET I M£ 

PRINT  71, VAR 
PRINT  69 

PRINT  68, (L# 1=1. 20) 

PRINT  77, (HIST (L)*L- 1,20) 

RETURN 

66  FORM AT < F6 . 4 ) 

67  F0RMATt/l6X,41HCH!  SOtJARE  VALUE  FOR  THIS  DISTRIBUTION  IS.F10.4/5 

68  FORMAT (2014) 

69  F3RMAT(lfX,<6HHlSTOr;RAM  OF  NORMALIZED  DATa*  7=0  AT  H  ON  S C A l E / ) 

70  F0RMAT(2fX,17H  =  X?ECTtfj  VALUE  IS.F1C.4/) 

71  F3^MAT(26X.11HVARIA\'CE  I  S,  6X .  FlC  .  4/ ) 

77  FORMAT ( 20F4 , 0  ) 

998  CONTINUE 

•  END 
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Figure  12  an  example  of  the  coroutine  ou'.put,  nia  Logi-au. 

of  data  input  representing  the  normal  probability  density  function,  and 
histogram  of  the  normalized  data. 

III.  METHOD  OF  SIMULATION 


This  section  will  pertain  primarily  to  the  logic  and  procedures 
utilized  in  simulating  the  network  model.  Input  data  requirements  and 
description  of  output  is  deferred  to  the  following  section,  Input-Output. 
The  three  parts  of  the  simulation  model,  or  program  will  be  discussed 
with  reference  to  the  subprograms  previously  mentioned.  The  flow  diagrams 
are  given  as  Figures  13  through  15  and  the  computer  listing  of  the  main 
program  is  found  on  pages  45  through  51 • 

All  arrays,  variables,  and  sums  are  initialized  at  the  start  of 
the  program.  The  random  number  seed,  RN,  and  number  of  simulation  trials 
NSIM,  are  defined  by  the  first  network  data  card  read  into  the  program. 
Should  RN  be  zero  the  program  is  terminated. 

Activities  specifying  the  network  model  and  their  attributes  are 
read  into  the  program  and  stored  in  the  array  ST0RE(N,  I).  Attributes 
of  each  activity  are  I-J  node  designations,  codes  for  the  probability 
density  functions,  parameters  of  the  density  function,  and  probability  f 
traversal  given  that  the  I  node  has  been  realized.  Subroutine  BETAD  is 
utilized  to  transform  the  PERT  a,  m,  and  b  time  estimates  into  parameters 
oi  and  Yof  the  Beta  probability  function  when  required.  As  the  activities 
are  read  into  the  program  the  number  of  activities  emanating  from  a  speci¬ 
fic  node,  NQ(I),  is  calculated. 


Nodes  and  their  logic  related  code  are  defined  on  the  next  inout 

the  variable  M”Cr2(J).  The  node  information 


cards  and  stored  as  values  of 
is  then  printed  as  output. 


Statement  79 


Sample  f 
distribution 


Probabilistic 


AND 

Deterministic 


Dump  emanating  NQ(j] 
activities  into  SET* 
schedule  each  activ- 
ity  * 

t±i  TNOtf  +  DIST 


XX=PAND(RN) 


Sim  attribute 
probability  of| 
acceptance 


schedule  one 
activity  from 
node  t  put  in  S5T$ 
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DATE  03/39/66  at  10582J  PAGE  NO. 

PROGRAM  r*AC 
COMMON  Rf 

COMMON  STATSflrOn#4) 

COMMON  StTS(  <0  »5 )  *  I  \f  its#  OUT # ^RA-MK { 1 )  #  inmc  d  ,  vi_c(  1 )  •  mf  a 
DlMtFS !0V  Tc  t 4 n 0 ) » NFL AG ( 4 ) , N«ST ( * ) 


DpSfcNSI 0*  A(  *  > .  X  C 40  0 .  A)  ,  'M(  100 )  ,  JflC l?Q ) 
01  MeNS I Ov  \3?’E0(  100 )  * mnQOE  <  3.00 ) 
DlVrMSIOV  NCRi  VClOO)  *STOREC^?IO»9) 


1 


INITIALISE  all  ARRAYS  amp  Sums 


1  CONTINUE 
PRINT  114 
PRINT  120 
RNr.07? 

CALL  RA:MU(RN) 
VT I ME=  0.0 
ETlMEsO.n 
JMAXsO 
M=  p 

DO  296  1=1,400 
00  296  J=l,« 
296  STW.<  l ,  J )  O  .  3 
DO  299  1=1,60 
TEC ITs-1.0 
JQ(  l )  =  0 
noc  j  >=n 
MNQOF  (  U  =  0 
NCRITC1  )  =  0 
DO  299  J=l,*i 
299  X<!,J)*0,0 


READIES  ! \  DATA  IN  l v ! T I AL 1 7 AT  10  ! 

RN  IS  RAVOO-1  RJMRER  SPED  A^D  NSIM  IS  MUMPER  OF  T&  1  AL  SIMIJL&T 

RF  AD  101,  RN,  'ISM 

IF{fm«e*,R68,T 

3  PRINT  101  ,  R  n,  NS  I v; 

C  PJT  DATA  O'*  TA^c 

DO  50  “si  ,  4 C  n 

READ  126, (A( J ), 1=1,7) 

IF C AC  1) -99. ) 4R, 51, 31 
49  NsMM 

A  (  p )  =  4  (  7  ) 

L  =  A  ( 5  ' 

GO  r0{1j|45i  1  ,  ^ -1  »  ^5*  i5  ,  I  1  L 

4  SlSsSO  -'T^C  A  (4)  > 

AC4)sA(3>-l.7«?S?l*5lG 
A  ( ? )  =  A  (  3  )  ♦  1 . 7  5  2  P  5  ?  *  <  1 1? 

I  r  (  A  t  4  )  )  0  5  1  ,  *o  ,  “  S 

951  £(0:0 
GO  T  O 

4  4  R  s  (  A  (  4  )  -  *.  (  3  )  1  /  (  A  (  b  )  -  '  (  U  ) 

CAI.L  iHTADC  O  A«_34A,  ) 

A ( 7 ) =  \ (  .0 


ooo  r»  cs  o  oonn 


FTN  1.4 


IT  105323 


Date  03/30/66 

AC6)=A( 3> 

AC3)*ALPHA 
A( 4 ) sGI HMA 
60  TO  46 
45  CONTINUE 
AC6)*1. 

A(7)sl. 

'  46  N3<A(1>)sN0ia<1>>*1 
DO  555  1  =  1.8 
555  ST0R0(‘M)  =  A(I) 

J»UX=M4X1F{ JMax, A<2)> 

50  C0MriN'J5 

51  co*Tms 
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PAGE  NO.  2 


READ  IN  NODc  INFORMATION 

COOE  FOR  MODES  EXCL-GR  :)ETs  l  ,  EXCl*0R  PR93=  ?  •  INCL-OR  DET 
INCL-OR  PR03=  4  #  A-iO-OFTs  5  ,  AND-PR0B=  6  . 


PRINT  130 
PRINT  131 
DO  65  m*1,M 

Pc AD  l?7,<NCRiT<!I)»l|*l»P4) 
1FCNCRIT(1)-9^)o?.66,6? 

62  PRINT  123, ( 4CRlT(i),?sl,?3»2> 
PRINT  12°, <\CR! T( I ) , I s?»?4»2> 
PRINT  131 

DO  63  Ml! si ,?4,2 

63  MNr>DF<NC«tT(«I  I )  T<  NJ I  ♦  3L  > 

65  CONTINUE 

IF (N-3? )66, 66, f 1 
61  PRINT  114 
PRINT  1.31 

66  no  67  Lsl  ,6'J 

67  NC«IT<L>s0 
DO  2*0  Lsl, 4 
NFlAGCUsO 

280  N^ST  <  L )  s  1 


NFLAG  DENOTES  TME  N0Uc  being  analysed,  N A X I m J f 4  OF  4  NODES 
READ  10i>,  (\’FLA3( !  > .  1  =  1 ,4) 

I  MI T I AL I ?E  SLTS 


J5  =  1 
INM( 

I N !  T  9  s  1 
CALL  SET ( JS ) 
NXsl 
JS  =  1 

KR\NM  J5>  =  1 
TAKE  Dir:.  ~-i'  A  r  Ar'  = 
DD  6“'  '  =  :  ,  V. 

n:i  6  ■  0  <  =  : ,  - 
60  0  A {  K  )  =  ST  0r,T  (  I »  <  ) 

MAAs A  {  1  ) 

N  C  T  :J  ’  ■  f  s*> 


C 


I  N  j  '0  =  -f  f.  <  I 


AT  105323 
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FTN  1. 

I 


3 


I 


4 


DATE  P3/3n/66 
IFCNAA-mOO. 53,52 


52  MsMAA-1 

DO  5B  K=l.NA 

58  NC0JNT  =  NrO*JMT*vi3(sO 

59  KArNCOU 4T ♦! 

53  ircXC\A»l>)ino,55,55 

55  NAsMA^l 
GO  TO  53 

56  DO  57  L=i.B 

57  XCNA.L)=A(L> 

60  COM  I  .«•*£ 

CHECK  RjrFER  AMD  NO(I>  F-< 03  NAIM 

PRINT  103.  C  <X(M,iO,K  =  l#  =  1  #M> 


PRINT  106.  M 
IF(M-15)015,915.916 
916  PRINT  114 
PRINT  141 

915  CONTI N-Jc 
C 

C  START  TRIAL  SIMULATION 
DO  995  KKrl.NS!  I 
C  INITIALIZE  TP  l  At.  SIMULATION 
DO  991  1=1.60 
T6(I)s-l.O 
991  J3( I ) -0 

DO  75  Ksl.M 
L  =  X('<.2) 

LLl*HNnDP(L) 

GO  TO  ( 11, 11.11. 1.1»  17.12 >.ULL 

11  JQ(L)=1 
GO  TO  75 

12  J3<L)  =  j:ja>*l 
75  COM  IN  1  = 

KRANK ( 1 ) =3 
TNOw=0, 

N00c= 1 

m 

j=i  1 

C 

C  RELEASE  *M30P»  SCHEDULE  4Ll  -MAM* TING 

C 

79  XVOOEs-'OI'E 
DO  155  t*l.  -i 
C 

C  DJMJ  SMAf'ATp.R  XCT  I V 1  T  I  I  ill)  SMS 

c 

15?  I“(X(L.U*X  ■<0D=)t55.5,165 

5  NNsNM  ’’D'lb)  4,|-“1 
LlL  =  mn (  J ) 

GO  in  ( 153, 4, 16 j, 6. 1 63, 6> ,LLl 

6  XXsRV/M 

S  J '  *  =  *1 . 

HO  -3  4  =  L  .  v 
SU”  =  ?M*  '♦%(<.  6  ) 

7  Sr  T  S  (  1FA.1  )  =  y  (  <  .  I  ) 


» 

activities 


uou 


FTN  1. 
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4  DATE  03/30/66  A7  105423  PAGE  NO.  4 

S6rSCri>'A,2)sX(<»?) 

S6T5<flFA,3)aTN3**D13T<X(<,3)»XC<.4>#X<K#5),X<<#6>#XCK.7)> 
CALL  SET(JS) 

GO  TO  93 
8  CONTINUE 

153  CONTINUE 
00  154  <=L#N.N 
SeTS<'IF  h#1>  =  X<<#1) 

St T3f  Hr  A#?)s *<(<#?) 

SETStlFAtSJsT-OA^OISTtXtK.SJrXtK.AJ.XtK.Sj.XtK.ftWXCK.?)) 
CALL  SET  ( J3) 

154  CONTINUE 
GO  TO  90 

155  CONTINUE  .  . 

UPDATE  CLOCK  AMU  09A4  NEXT  EVENT 

90  TM04=SSTS(MLC<1>»3> 

JaS5TS(“.LCil).2) 

IsStTSMLC(i).l) 

OUTM  • 

CALL  SET C JS ) 

IFtJ-I  >95,100.1.9? 

192  !F( JO<  J) ~1) 91 , 92, 93 

95  00  99  un#,j 
NlJsX(IJ#l) 

IF<N!J-J)91,97,94 

96  IF(NIJ-I)97,94,91 
94  N39IT( l J)aNC9IT( T J)+l 

97  LlLsNNOOECMIJ) 

GO  TO  (15#15»15#15#l6#16)  #LLl- 

15  J?<NIJ)«1 
GO  TO  94 

16  J0( MI J> * J0( 4 1 J > ♦! 

9fl  CONTI  Nile 

GO  TO  91 
9?  N9f?cO(J)sI 

91  J0(J)=n 
N00c=J 
T5(J)aTV0W 
IF(N0(  JJU(l'),i?0. 79 

93  JOf  J) s J  )( J) *1 
GO  rn  no 
101  P^INT  199 
199  F0f?NAT(5lir.9  n  O 
P9INT  106,  l, J 
GO  TO  1 
320  CONTINUE 

ET!4r.st:nM  =  -rc(  J) 

V  T I  4E  =  '#TIK=»rs(  J)*Tr;(  J) 

no  22?  1^=1#  » 

lr*(  VH  i!  H)  >2^2.  A?  >,?  M 
2?  s  i  “  '  7 ~  ( *  -  L  AT,  {  l  -i )  )  >  >P  > ,  ?  >  ) ,  ??  ) 

229  'T=I~*l 

N r» ?:  f  <  H)s\MST(  H)*l 
ST  A  re  (  V*GT(  H)  ,  H)sTE(  NFL  AG<  N) ) 


FT 
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222  CONTINUE 

ESTABLISH  CRITICAL  EVENT 
% 

MJUsJ 
NNsN 

201  Ls'Jp-POC:^) 

.  XL=L 

DO  215  K=1.NN 
NRORsNN-K+l 

IF(X(NPO,J,l)-XL)2C5,2i?,205 

202  IF  (X(NO0*,  2) 203.205 

203  NC»lT<*Rnw)sMCRirCNROw)*l 

205  CONTINUE 
KNsL 

IF<L-li2r.6. 206,211 

206  CONTINUE 
IF<«LC> JS>>995,995,905 

905  INITSsJ 

CALL  SET  <JS> 

995  CONTINUE 

CALCULATE  STATISTICS 

XMSJMsNSfM 
ETIMFs  =  T»H5/<i*SlM 
VT I  Mr*V  r  I  ME/  <  NS  Il-rfl  *16>ET  IrtE 
PRINT  10R.5THS 

PRINT  10O,VTIM? 

C  CRITICAL! IV  CALCULATIONS 

PRINT  111 
DO  200  1=1 .  J 
CR  t  7*  NCR l T ( I ) 

RsCRI  T/X'iSri 

PRI.NT  111,  XCI,  1),  X(1, 5).  *3 

2on  continue 

PRINT  H4 
NNsO 
C 

C  HISTOGRAM  CAlC'JLAT  i  0  mS 

00  422  !Msi,t 
IF(NFLA3( IH))422, 422*423 

420  NT* l H 
PR  I  iT  131 

PRINT  U2»  N"L AC(  I H) 

PRINT  t3l 

PRINT  1  j>,  MIST  (  H) 

XHsiPSTUH) 

TTsXH/VNSlv 
PRI’jT  1  3  < ,  T  T 
Call  HI  iT(vr,  4  J  S  T  <  I-i)  ) 

I*  (I-*-?)  i??,  «2l 

421  PRINT  114 

422  COUINJs 
GO  n  1 


■yfwp*  *» _»>  *  .  ■  ,  •  -,  ,  ■ 


'* 


FT#  1.4 


« 


50 

DATE  03/30/66  AT  105323  RASE  NO.  6 

887  CONTINUE 

888  CALL  E/IT 

101  FORMATIFS.J, !«> 

102  FD=MAT(4F3.a,F5.0> 

103  F09M4  f < 5F 10. 2) 

106  FORMAT 121 14) 

108  F0»MAT(2SX.27H=<:,ECT  =  D  T I  HE  LAST  EVENT  |S,F10.4/> 

109  F  QP.M1  T  ( 26X,  2  7H  VARIANCE  TIME  LAST  EVENT  IS.F13.4/) 

110  FOR-iAT  ( /31X , 1 7HCRI T I  CAL  I  TV  INDEX. /3J  X, 4H  I,4H  J.10H 
1) 

111  FORMAT<33X,2F4.0,F10.2> 

112  F0RMArf27x,2lH*HlST03RA1  DATA,  NODE. 14, IN*. /) 

114  FORMATU*l) 

115  FORMA  T</20H> 

120  FORMAT  </27<*25H*3ERT  NETWORK  SIMULATIONS//) 

126  FORMAT ( 2F5. •). 2?10 •  ?. 3F5. 3 ) 

127  F0R1ATI24I3) 

128  FORMAT C 15X, 5M*M0-1E. 2X,  1213, lw* ) 

129  F0Q.MAT(18X,3R*TT»E,2X.12I3»1H*) 

130  F0RMAT<34X,1?H*IV?UT  »)ATA«> 

131  F0*MAT(//> 

132  FORMA T ( 2*3 X » 27MN\JM35ri  OF  MOD :  REALIZATIONS,  1 10  ) 

133  FORMAT(/,2n.<,29HaR0dARlLlTr  OF  REALIZATION  IS.F10.4/) 

999  F0RMaT(5F10.4) 

End 


CR 


FT*  1,4 


DATE  P3/30/66  AT 

FUMCMOV  3AmD(X) 
X*X»1.0=*8*23 
I-X*l.PE-6 

y*i 

xs(x-r-v*i.0P4.fu*i.9E-fi 

DA'JDsX 

RcTJkM 

EM3 


10W3 


CRJT/ 
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A  maximum  of  four  nodes  may  be  flagged  for  analysis  during  the 
simulation  of  the  network  by  including  them  on  the  next  data  input  card. 
These  node  numbers  are  retained  as  values  of  the  variable  HFLAG(XH). 

A  buffer  array,  X(M,  K) ,  is  utilised  as  the  permanent  storage  area 
of  the  activities  and  their  attributes,  which  are  transformed  and  ordered 
by  1  node  designations  from  the  array  STORE  (N,  I).  The  contents  of 
X(K,  K)  ate  then  printed  as  output. 

The  basic  procedure  in  simulating  the  network  is  one  of  working 
through  the  network  from  node  to  node,  releasing  activities  and  schedu¬ 
ling  them  in  the  master  time  file  SETS(J,  K),  and  recalling  the  activities 
from  SETS(J,  K)  as  they  have  been  completed. 

The  actual  simulation  of  a  network  can  be  more  easily  explained 
by  tracing  through  a  simulation  cycle. 

Each  cycle  begins  with  initialisation  of  TE(I),  time  node  I  is 
released;  TNOW,  current  time  of  system;  and  JQ(J),  required  number  of 
activities  terminating  at  a  node  necessary  for  realization  of  the  node. 
JQ(J),  for  all  nodes,  is  set  equal  to  1  for  Exelusive-or  and  Inclusive- 
or  type  nodes,  and^  1  for  AND  type  nodes.  The  start  node  is  then  real¬ 
ised  by  setting  NODE  =  1,  the  first  node  of  the  network. 

As  each  NODE  is  realized,  control  is  transferred  to  statement 
79t  which  causes  the  emanating  activities  to  be  released  and  scheduled. 
Should  the  output  node  types  be  deterministic,  all  activities  emanating 
from  that  node  are  scheduled  and  placed  in  the  master  time  file  SETS (I,  K) 
with  scheduled  time  equal  to  i’NOW  +  DIST,  As  previously  mentioned, 
the  variable  DIST  takes  on  a  time  value  as  a  random  variable  of  the  prob¬ 
ability  function  for  each  activity  released. 
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Should  the  output  node  type  be  probabilistic,  a  random  number  is 
generated  and  used  as  a  reference  in  determining  which  activity  of  the 
activities  emanating  from  the  node  should  be  scheduled.  Only  one  activ- 
.  ity  is  scheduled  with  time  equal  to  TNOW+  DIST;  the  others  are  dropped. 

Each  act  of  placing  an  activity  in  SETS  is  accompanied  by  a  call 
of  Subroutine  SET(l) ,  which  causes  ordering  of  the  activities  in  sets 
through  a  marker  arrangement  23  previously  mentioned. 

After  the  appropriate  number  of  activities  have  been  released  and 
filed  in  SETS  (I,  K),  control  is  advanced  to  statement  90,  which  causes 
the  next  event  to  occur.  The  activity  filed  in  SETS  (I,  K)  marked  for 
next  exit  by  HLC(l)  is  drawn  from  SETS  (I,  K)  and  the  system  time,  TNOW, 
is  updated.  TNOW  is  set  equal  to  the  scheduled  time  of  completion  of  that 
Activity. 

A  comparison  of  the  I  and  J  attributes  of  the  completed  activity 
is  made  to  see  if  a  feedback  condition  exists,  where  J <  I.  If  not,  then 
JQ(J),  number  of  completed  activities  required  node  J,  is  investigated. 

If  JQ(J)  >  1,  then  JQ(J)  =  JQ( J)-l,  which  means  more  activities  have  yet 
to  be  completed  before  node  J  can  be  realized.  Control  is  returned  to 

1 

statement  90  for  next  release  of  a  scheduled  activity  from  SETS(l,  K), 

If  JQ(J)  -  1,  node  J  is  released  and  the  time  of  release,  TE(J), 
is  set  to  TNOW,  Precedence  for  criticality  measurement  is  established 
by  setting  NPRED(J)  =  I. 

The  variable  NQ(J),  or  number  of  activities  emanating  from  node 


J,  is  investigated  to  determine  if  an  exit  node  has  been  realized.  If 
NQ(J)>  0,  control  is  transferred  to  statement  79  with  NODE  =  J  for 
release  of  node  J  as  described  above,  ’’.hen  J)  =  0,  the  network  has 


been  completed  since  J  was  a  terminating  node,  and  statistics  are 


\  — s  “  •,/■ 
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Control  is  returned  to  the  beginning  of  the  simulation  phase  for 
another  network  simulation.  This  process  continues  until  the  required 
number  of  simulation  trials,  NSIM,  has  been  attained.  Final  output  sta¬ 
tistics  are  then  computed  and  printed. 

Intentionally  emitted  from  the  preceding  discussion  is  the  special 
case  of  feedback,  where  J  K  I  for  a  completed  activity.  Feedback  returns 
network  control  to  a  node  previously  realized  in  the  network.  A  meaning¬ 
less  situation  could  exist  since  scheduled  activities  positioned  between 
nodes  I  and  J  may  not  have  been  completed  and  are  still  in  the  master 
time  file.  The  network  of  nodes  and  activities  comprising  a  "feedback 
loop  is  depicted  in  Figure  l6. 


FIGURE  16 


FEEDBACK  NETWORK 

Nodes  4,  5,  and  6  define  the  feedback  loop,  and  the  node  numbering 
is  such  that  J<  NIJ  ^  I,  NIJ  represents  the  nodes  positioned  within  the 
loop. 

Another  situation  that  could  exist  stems  from  the  fact  that  each 
JQ(J)  within  the  feedback  loop  has  been  set  to  zero  as  the  nodes  were 
released.  This  programming  aspect  is  corrected  by  re-establishing  JQ(J) 
for  each  node  NIJ  and  I  within  the  feedback  loop,  J<  NIJ  <  I.  Thus, 
recycling  is  permitted  for  this  feedback  portion  of  the  network. 

Caution  is  advised  in  designing  the  network  model  to  eliminate 
undesirable  feedback  situations.  For  example,  consider  the  network  shown 
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from  a  node  within  the  loop  to  the  external  network. 


As  the  network  is  constructed,  the  first  realization  of  node  6 
would  release  activities  6,  8  and  6,  10.  The  network  beyond  node  10 
would  then  be  activated.  Additional  realizations  of  node  10  by  repeated 
looping  among  nodes  4,  6,  8,  4  could  not  cause  another  release  of  node 
10  since  JQ(lO)  =  0  after  the  first  realization.  However,  should  the 
analyst  desire  to  halt  realization  of  the  external  network  beyond  node 
10  until  activity 8 , 14  is  released,  a  restraining  activity  14,  15  could 
be  constructed.  Node  10  can  have  only  one  terminating  activity  for  this 
scheme  to  work  properly. 

Consider  a  more  complex  network  than  that  discussed  above  where  a 
second  feedback  loop  is  compounding  the  effects  of  the  first.  The  network 
is  shown  in  Figure  18. 


TWO  LOOP  NETWORK 


Since  JQOO,  for  J  <  K  <  I  is  reset  lay  the  program  in  a  feedback 
loop,  .repeated  looping  from  nodes  4,  6,  8,  4  will  not  cause  more  than  one 
release  of  node  10  due  to  this  loop  alone.  However,  should  activity  10,  2 
be  released,  feedback  looping  would  occur  in  both  loops.  Restraining 
activity  15,  16,  inserted  by  the  analyst  would  prevent  the  possibility  of 
repeated  releases  from  node  8  to.^ihe  external  network,  since  each  real¬ 
isation  of  the  loop  10,  2,  4,  6,  10  would  cause  node  8  to  be  realized. 

IV.  INPUT-OUTPUT 

An  important  attribute  of  the  computer  program  is  minimization  of 
user-time  in  preparation  of  input  data  for  the  simulation  program.  This 
section  describes  the  input  data  cards  for  computer  processing  of  input 
data  and  program  output.  A  sample  listing  of  input  data  and  output  is 
provided  for  reference  on  pages  57  through  60. 

Only  one  input  card  per  activity  is  required  containing  I-J  node 
designations,  parameters  or  time-estimates  describing  the  probability 
density  function,  cede  for  the  density  function,  and  probability  of 
acceptance  should  the  I  node  be  realized. 

Codes  for  the  probability  density  functions  and  input  parameters 
or  time  estimates  are  listed  as  follows. 


Density  Function 

Code 

Parameters  or  Estimates 

Exponential 

'  1 

u 

Normal 

2 

2 

U  o 

Rectangular  Discrete 

3 

«  2 

U  W 

Constant 

4 

T 

Beta 


5 


a  m  b 
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*GERT  fcEUORK  SIMULATION* 


*4771000 


♦INPUT  DATA* 


♦NODE 

1  2 

3 

4  5  6 

7  -0  -0 

-0  -0  -0* 

♦TYPE 

5  2 

2 

5  2  2 

1 

2  -0  -0 

-0  -0  -0* 

# 

1,00 

2.00 

-0.00 

-C.00 

4.00 

1.00 

1.00 

1.00 

2,00 

3.00 

2.00 

-  G ,  0  0 

4.00 

1.00 

1.00 

0.2o 

2,00 

5.00 

2.00 

-0.00 

4.00 

1.00 

1.00 

0.50 

2,00 

7.00 

2.00 

-0.00 

4.00 

1.00 

1.00 

0.3Q 

3,00 

4. no 

-0.00 

-0.00 

4.00 

1,00 

1.00 

Q.lQ 

3.00 

5,00 

2.00 

-0,00 

4.00 

1.00 

1.00 

0.30 

3,00 

6.00 

-0.00 

-0,00 

4.00 

1.00 

1.00 

0.30 

3,00 

7.00 

2.00 

-c.co 

4.00 

1.00 

1.00 

0 . 3  0 

4,00 

3.  GO 

1.00 

-0,00 

4.00 

1.00 

1.00 

1.00 

EX°EC  f  £.'»  Thf  LAST  EVENT  is  2,293" 
VARIANCE  TH:  LAST  EVENT  IS  0.597? 


CRITICALITY  I  \'UEX 

1  J  CriU 

1  2  1.00' 

?  3  0.21 

2  5  0.48 

2  7  0.31 

•3  4  0 , 0  0 

3  5  0.08 

3  6  0.08 

3  7  0.05 

4  3  0.04 


i 

i 
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♦histogram  data,  node  4* 


NUMBER  of  MODE  REALIZATIONS  29 

PROBABILITY  OF  REALIZATION  IS  0.0290 


EXPECTED 

TIME 

OF 

COMPLETION  this 

node 

IS 

2.2759 

0 

0 

variance 

time 

OF 

completion  this 

node 

IS 

0,4067 

g  2 

3  4  5 

6  7  9 

9 

10 

11  12  13  14 

15 

16 

17  18  19 

20 

S23 

5  0  1 

0  0  0 

0. 

0 

0  0  0  0 

0 

0 

0  0  0 

0 

o  big  small  range 

o  5,000  2,000  3.000 

0 


♦HISTOGRAM  data,  node  5* 


NUMBER  or  NODE  REAL  I /AT  IONS  563 

PROBABILITY  OF  REALIZATION  IS  0 . 5*30 


2  3  4 

60  0  69 

BIG 
6,000 


EXPECTED  TIME  OF  COMPLETION  THIS  NODE  IS  '  2.3321 

VARIANCE  TIME  OF  COMPLETION  THIS  NODE  IS  0,7014 

5  6  7  3  9  10  11  12  13  14  15  16  17  18  19  20 

940109  ,  000  0  000000 

Small  range 

2,000  6,000 


I 
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•HISTOGRAM  DATA.  NODE  6* 

NUMBER  OF  NOOt  REALIZATIONS  80 

PROBADILITY.OF  REALIZATION  IS  0.0600 

EXPECTED  TIME  OF  COMPLETION  THIS  NODE  IS  2.1125 

VARIANCE  TIME  OF  COMPLETION  THIS  NODE  IS  0.0995 

2  3  4  5  6  7  3  9  10  11  12  13  14  15  16  17  13  19  20  2 

71  90000000000  0000000 

816  small  ranss 

3,000  2.0PO  1.000 

•HISTOGRAM  OATA.  NODE  7* 

\ 

NUMBER  0“  NOOE  REALIZATIONS  357 

PROBABILITY  OF  REALIZATION  IS  0.3570 

EXPECTED  TIME  OF  COMPLETION  THIS  NODE  IS  2,2717 

VARIANCE  TIME  OF  COMPLETION  THIS  NODE  IS  0.5343 

5  6  7  3  9  10  11  12  13  14  15  16  17  18  19  20 

5100000000000000 
SMALL  range 

2,000  4.000 


2  3  4 

312  0  39 

6IG 
6,003 


5 
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Tina  estimates  a,  a,  and  b  for  the  Beta  Distribution  are  converted 
at  input  time  through  SUBROUTINE  BET  AD  into  parameters  and  of  the 
Beta  probability  density  function. 

The  data  input  cards  for  each  activity  are  punched  in  FORMAT 
(215,  2F10.2,  3F5.0),  in  the  following  sequence  for  distribution  codes 
1-4, 

I  J  u  <T  Code 

and  in  the  following  sequence  for  distribution  code  5 

I  J  a  a  Code  b 

GERT  logic  node  types  and  codes  are  given  below. 


INPUT  OUTPUT  CODE 

Inclusive -or  deterministic  1 

probabilistic  2 

Exclusive -or  deterministic  3 

probabilistic  4 

AND  deterministic  5 

probabilistic  6 


Each  network  node  from  node  1  to  the  end  nodes  is  defined  with 
its  respective  node  type  by  input  to  the  program  in  FORMAT  (2413)  as 
indicated  below. 

Node  Type  Node  Type  Nod6  Type 

1  5  2  5  3  4 


As  many  data  cards  are  utilized  as  necessary  to  input  the  nodes 
and  types. 

Nodes  that  are  to  be  analyzed,  causing  a  complete  print  out  of 
statistics  gathered,  are  denoted  as  FFLAC*  ar.i  i.-:  FC 


:  > 
'  / 


The  data  deck,  along  with  the  necessary  control  cards  are  listed 


e 
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below: 

Card  Column 

12345678  9  10  11  12 
004771000 

I  J 


Random  number  seed  and  NS1M 
Activity  cards,  one  per  activity 


0  9  9 

001005002  0  0  5 

0  9  9 

00060007 


JL 

End  of  activity  list 

INode  numbers  and  types,  12  sets 
per  data  card 

End  of  node  list 

Nodes  to  be  analyzed 


As  many  networks  may  be  analyzed  as  desired  in  one  computer  run  by 
providing  the  sequence  of  data  cards  described  above  for  each  network. 

0 

The  network  sets  of  data  cards  are  stacked  one  on  the  other  with  a  blank 
card  placed  at  the  end  of  the  last  network  to  signify  end  of  program  data. 

Should  the  Chi  Sjquare  Goodness  of  Fit  Test  be  required;  ten  input 
cards,  each  containing  the  expected  cell  block  size  obtained  from  the 
tabulated  probability  distribution  function  utilized,  is  entered  in 
FORMAT  (F6.4).  The  computer  program  is  designed  for  symmetrical  distri¬ 
butions,  so  the  ten  equally  spaced  cells  would  represent  the  area  under 
the  density  function  from  the  mean  to  one  tail. 


Computer  output  is  divided  into  five  sections:  (l)  input  data, 

(2)  expected  and  variance  of  network  c caplet*  ;r.  tire,  (3)  criticality 


index,  (4)  node  statistics,  and  if  required  (5)  chi  square  test  results. 
Discussion  follows  based  on  this  division. 

in  echo  check  of  the  random  number  seed,  and  number  of  simulation 
trials  precedes  a  complete  print  out  of  the  activities  and  their  respec¬ 
tive  attributes.  Each  activity  contained  in  the  buffer  storage  array 
is  printed  in  the  form  utilized  during  simulation,  as  follows  for  dis¬ 
tribution  codes  1-4. 


I.  J  u  code  -  -  p^ 

For  activities  distribution  code  5*  the  following  is  printed  out. 

I  J  **  &  code  a  b  p^ 

The  next  section  provides  the  expected  and  variance  of  network 
completion  time  regardless  of  the  number  of  end  nodes. 

A  criticality  index  follows  for  all  activities  in  the  network  as 
the  relative  frequency  of  occurrence  of  an  activity  on  the  critical  path 
of  the  network.  Activities  possessing  the  I  >  J  feedback  characteristic 
•  will  have  recorded,  not  criticality,  but  total  number  of  times  realized 
divided  by  number  of  simulation  trials  since  many  utilizations  of  a  loop 
may  occur  during  each  simulation  trial. 

For  each  node  flagged  at  input  for  analysis,  the  output  will  con¬ 
sist  of,  number  of  realizations,  probability  of  realization,  mean  and 

variance  of  realization  time,  histogram  of  these  times,  high  tine,  low 

* 

time,  and  range. 

The  Chi  Square  Goodness  of  Fit  Test  is  printed  out  only  when 
Subroutine  CHI  is  utilized  and  theoretical  distribution  input  provided. 
The  calculated  chi  square  value,  mean  and  variance  of  the  sample  dis¬ 


tribution,  and  histogram  of  the  nor-.alir,?l  distribution  i.v  provided. 


t*;*cv*  ftofnr  t *« 


CHAPTER  IV 


C 

(.. 


APPLICATIONS  CP  THE  SIMULATION  MODEL 

The  purpose  of  this  chapter  is  to  present  several  applications 
of  GERT  network  analysis  through  use  of  the  GERT  Simulation  Program. 
Several  example  problems  have  been  chosen  to  demonstrate  general  appli¬ 
cability  of  the  computer  program  to  various  type  networks. 

Three  sample  probabilistic  networks  are  simulated  and  the  results 
compared  with  analytical  solutions  such  as  from  the  GERT  computer  program 
described  by  Pritsker.^  The  aim  is  not  to  perform  a  detailed  analysis  of 
the  individual  networks;  moreover,  to  indicate  the  approach  to  problem 
solution  with  the  GERT  Simulation  Program  and  to  provide  a  measure  of  the 
validity  of  simulated  results. 

One  other  network,  consisting  of  AND  nodes,  is  considered  to  in¬ 
vestigate  assumptions  developed  in  the  literature  on  PERT  networks.  The 
resultant  distributions  are  analyzed  using  normally  distributed  and  beta 
distributed  activity  times.  The  Chi  Square  Goodness  of  Fit  Test  is 
applied  to  test  the  hypothesis  of  normally  distributed  network  realiz¬ 
ation  times. 

I.  PROBABILISTIC  NETWORKS 

The  first  network  depicts  an  analysis  of  a  Research  and  Develop- 

■^A.  Alan  B.  Pritsker,  GERT:  Graphical  Evaluation  and  Review 
Technique ,  The  RAIID  Corporation,  RM-h973-NASA  (Santa  Monica,  April,  1966), 
Appendix  A. 
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sent  project  given  In  a  recent  article  by  Graham,  and  discussed  in  rela¬ 
tion  to  GERT,  by  Fritsker.^  The  network  is  constructed  in  Figure  19 
For  each  branch  or  activity  of  the  network,  the  probability  that  the 
branch  is  realized  given  the  preceding  node  is  realized  is  depicted  by 
p^,  and  the  constant  time  associated  with  activity  k  is  depicted  by  t^. 
Bote  that  there  are  two  terminating  nodes  of  this  network;  success,  s; 
and  failure,  f. 

The  criticality  index  of  the  relative  frequency  an  activity  con¬ 
tributed  to  network  realization,  is  presented  in  Table  I.  The  ability 
of  the  simulation  model  to  select  activities  emanating  from  a  probabil¬ 
istic  output  node  can  be  checked  by  comparing  the  probabilities  associated 
with  activities  1-2  and  1-6,  also  activities  6-9  and  6-8,  in  Figure  19 
with  the  criticality  index  of  the  activities  in  Table  I.  Activities  1-2 
and  1-6  were  defined  with  probabilities  of  0.?  and  0.3  which  compare 
favorable  with  simulated  criticality  of  0.68  and  0.32,  respectively. 
Similarily,  activities  6-8  and  6-9  defined  with  probabilities  of  (.3) (.5) 
each  compare  with  simulated  criticality  of  0.1?  and  0.l6.  The  slight 
difference  in  each  case  can  easily  be  attributed  to  sampling  error. 


rearson  Graham,  "Profit  Probability  Analysis  of  Research  and 
Development  Expenditures",  The  Journal  of  Industrial  Engineering.  Vol. 
XVI,  No.  3*  May  June  1965»  pp.  I86-I9I. 

^Pritsker,  Oo,  Cit. ,  p„  57. 
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TABLE  I 


SIMULATED  CRITICAIITY  INDEX  ON  ACTIVITIES  IN 
THE  RESEARCH  AND  DEVELOPMENT  PRC® LEM 


Activity  I 

J 

CRITICALITY 

1 

2 

0.68 

1 

6 

0.32 

2 

3 

0.18 

2 

4 

0.15 

2 

5 

0.24 

2 

7 

0.10 

3 

12 

C.18 

4 

12 

0.15 

5 

12 

0.24 

6 

8 

0.17 

6 

9 

0.16 

7 

17 

0.10 

8 

10 

C,09 

8 

11 

0.07 

9 

17 

0,l6 

10 

17 

0.09 

11  . 

16 

0.07 

12 

13 

0.57 

.  13 

14 

0.40 

13 

15 

0.17 

14 

16 

0.40 

15 

17 

0.17 

Node  analysis  was  conducted  on  trades  l6  and  17  for  probability  of 
realization  and  resultant  distribution  of  realization  times.  A  test  of 
proportions  using  the  normal  approximation  is  utilized  to  test  the  hy¬ 
pothesis  that  the  simulated  probability  of  realization  is  equal  to  the 
analytical  probability,  at  the  5  per  cent  level  of  significance,  for 
each  possible  time  of  realization. 

The  test  statistic  is,  Z  =  (X  -  Np)/  /  Npq,  where  Z  is  the 
standardized  variable,  X  is  the  simulated  frequency  of  successes,  p  is 
the  analytic  proportion  of  successes,  q  is  the  proportion  of  failures, 
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q=l-p,  and  K  is  the  sample  size.  A  two  sided  test  is  applied  with  a 
critical  region,  -I.96  <  Z  <  I.96.  Results  of  the  node  analysis  are  as 
depicted  in  Table  II  where  the  Z  score  is  given  for  each  pair  of  simulated 
and  analytical  frequencies  at  the  discrete  time  intervals. 


TABLE  II 

RESULTS  OF  NODE  ANALYSIS  ON 
RESEARCH  AND  DEVELOPMENT  PR0BU5M 


Node 

Time  of 

Simulation 

Analytical 

Z  Statistic 

Realization 

Frequency 

Frequency 

11.0 

74.0 

60.0 

.904 

16 

16.0 

398.0 

411.6 

.  1 

$ 

Vn 

6.0 

158.0 

150.0 

.710 

* 

10.0 

104.0 

112.0 

vr> 

0 

GO 

• 

1 

17 

11.0 

91.0 

90.0 

.111 

16.0 

175.0 

176.0 

.116 

The  test  of  proportions  utilised  indicates  there  is  no  reason  to 
reject  the  hypothesis  that  the  simulated  and  analytic  results  are  equal, 
as  all  Z  scores  are  well  within  the  critical  region. 

The  second  network  considered  for  simulation  is  the  Thief  of 
Bagdad  abstracted  from  Parzen  and  analyzed  in  GERT  .  The  problem 
concerns  a  thief  in  a  dungeon  faced  with  the  selection  of  three  doors 


amanuel  Parzen,  Stochastic  Processes.  (San  Franciscos  Holden- 
Day,  Inc.,  19&2),  p.  50. 

^Pritsker,  op.  cjlt. ,  p.  45, 
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leading  to;  (l)  freedom;  (2)  a  long  tunnel  returning  to  the  dungeon; 
and  (3)  a  short  tunnel  returning  to  the  dungeon.  Experience  has  no 
effect  on  selection  given  he  makes  an  incorrect  choice  and  returns  to 
the  dungeon.  The  network  is  depicted  in  Figure  20  with  constant  activity 
time  assignments,  t,  and  probability  of  selection,  p. 


This  network  was  chosen  because  the  feedback  loops,  one  par/allel  to 
the  other,  would  provide  a  challenge  to  the  computer  program. 

The  expected  value  of  node  5*  for  the  1000  simulated  trials 
was  3*057  with  variance  equal  to  17.684.  The  criticality  index  of  .76 
for  activity  3-2  and  .77  for  activity  4-2  indicates  equal  probability  of 
selection.  The  histogram  of  realization  times  for  node  5  is  shown  in 
Figure  21, 

The  expected  time  to  realize  node  5  can  be  determined  by,^ 

E(t^)  =  tl2  +  +  l/p25(  P24t24  +  p23l23  ^ 


0  +  0+  .4(  .9+  .3  )  =  3*0 
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A  fnorable  comparison  is  noted  between  analytical  and  Simla  ted 
▼aloes  for  expected  tine  of  realisation,  node  5*  of  3*0  and  3-057* 

A  test  of  proportions  is  utilised  to  farther  check  the  accuracy 
of  the  s inula tion  results  for  probability  of  realising  node  5  in  less 
than  c  tifce  units  is 

P  s  p(t^<2)  =  p(t=0)  +  p(t=l)  -  A  +  (.4}(.3)  =  *52 

As  defined  for  the  previous  network,  the  test  statistic  is 
Z  *  (X  -  Ip)/  YHpq  •  Application  of  the  test  to  the  s  inula  ted  results 
of  • 513  analytical  results  of  .52  for  p(t^<  2)  at  the  5  per  cent 

level  of  significance  results  in  the  conclusion  of  no  difference  in  the 
two  procedures  for  this  tine  interval. 

The  bomber- interceptor  problem  is  chosen  as  the  third  example  of 

networks  containing  probabilistic  type  nodes.  An  interceptor  nay  kill 

or  br  lied  on  the  first  pass  at  the  bomber  or  it  nay  niss  and  nake 

another  pass  at  the  bomber,  again  facing  the  chance  of  killing,  be  killed, 

or  missing.  Successive  passes  are  allowed  with  the  chance  of  mission  abort 

after  each  p^ss.  The  network  is  depicted  in  Figure  22  with  probabilities, 

« 

p,  of  activity  realization  and  constant  times  associated  with  each  activity. 
Bote  that  a  small  feedback  loop  represents  the  successive  passes  of  the 
interceptor. 

The  criticality  index  for  each  activity  in  the  network  is  presen¬ 
ted  in  Table  III.  Comparison  with  the  probabilities  as  given  in  the 
network  indicate  only  slight  differences  attributed  to  sampling  error. 

The  criticality  index  for  activity  4-3  represents  total  number  of  real- 


FIGURE  22 

G2RT  NETV/CFK  FOP  BOMBER  FiTEPCEFTOR  PROBLEM 


lsations  over  tbs  1090  sinlition  trials,  and  is  used  as  a  s isolation 
diagnostic  9n d  not  a  saason  of  criticalness. 


TAME  m 


SH0LATED  CRmCAUTT 
ACTIVITIES  OF  TEE  BOMBER 


l  FOR 

*TCR  PROBLEM 


Activity  I  J 


1 

2 

2 

2 

3 

3 

3 

3 

4 


2 

3 
5 
7 

4 


7 

3 


.  CRITCCAUTT 

1.00 

0.21 

0.48 

0.31 

0.00 

0.06 

0.08 

0.05 

0.04 


Bode  Analysis  was  conducted  on  nodes  5»  6,  and  7,  for  probability 
of  realization  and  distribution  of  realization  tine.  Results  of  the  net¬ 
work  simulation  as  contrasted  with  those  obtained  through  application  of 
.the  analytic  G2RT  computer  program  are  summarized  in  Table  IV.  A  test  of 
proportions,  based  on  the  normal  approximation  and  Z  test  statistic  pre¬ 
viously  utilized,  is  applied  to  the  probability  of  realization  obtained 
for  each  node.  The  number  of  simulation  trials  was  1000. 
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TABUS  17 

RESULTS  of  uooe  analysis  ok  bomber  interceptor  problem 
FOR  PROBABILITY  AND  EXPECTED  TIKE  OF  REALIZATION 

Probability  of  Realization  Expected  Tine 


Bodes 

Simulated 

Analytical 

Simulated 

Analytical 

Z  Score 

5 

.56 3 

.567 

2.332 

2.2 48 

-0.19 

6 

.060 

.067 

2.112 

2.111 

1.64 

7 

.357 

.36? 

2.272 

2.384 

-0.65 

There  is  no  reason  to  reject  the  hypothesis  of  equality  of  simu¬ 
lated  and  analytical  probabilities  of  realization  at  the  5  per  cent 
confidence  level.  Although  the  Z  score  for  nod9  6  is  high,  it  falls 
within  the  acceptable  confidence  limits  oi*  Z  =  t  1.96. 

The  histograms  for  time  of  realization  of  nodes  5*  6,  and  7  are 
given  in  Figure  23  .  While  a  complete  cell  by  cell  comparison  of  expected 
versus  observed  histogram  cell  size,  utilizing  the  statistic,  would 
provide  the  strongest  tost  of  the.  accuracy  of  the  simulation  results, 
it  is  felt  that  comparison  on  tne  cell  from  each  histogram  would  provide 
sufficient  information.  The  test  of  proportions  is  applied  to  the  time  * 
of  realization  t  <  3  for  each  node  and  is  shown  in  Table  V . 
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TABLE  V 

ANALYSIS  OF  TIME  OF  R3AUZAXKW  (t<3) 

WXBS  4,  5,  AMD  7,  OF  THE  BOMBER  INTERCEPTOR  PROBISC 


Bode 

Simulated 

Analytical 

Z  Score 

5 

.480 

.50 

-1.26 

6 

.071 

•06 

1.46 

7 

.312 

.30 

.83 

Tbs  test  of  the  sample.  tine  values  for  all  three  nodes  indicates 
no  significance  at  the  5  per  cent  confidence  level.  However,  the  high 
score  for  nodes  5  end  6  indicate  nore  simulations  of  the  network  are 
necessary  to  establish  nore  accurate  results. 

II.  DETERMINISTIC  NETWORKS 


The  network  chosen  to  exhibit  solution  of  deterministic  GERT 
networks,  consisting  of  all  AND  nodes,  also  serves  a  dual  purpose  of 
investigating  PERT  networks  and  related  theory.  The  network  is  depicted 
in  Figure  24  with  activity  duration  times  given  as  the  PERT  a,  m,  and  b, 
tine  estimates.  Also  included  above  each  time  estimate  in  the  network 
diagram  is  the  mean  and  variance  calculated  from  the  PERT  approximation 
formulas: 

=  B(x)  =  (a  +  +  b)/6 

VAR(x)  =  (b  -  a)2/3 6  . 


This  network  was  simulated  using  each  of  two  probability  distri- 
butions,  Beta  and  Normal,  and  in  two  configurations  or  forms.  Form  1 
contains  the  network  and  activities  as  river,  by  Moder  and  Phillies  and 


1 


78 


7 

presented  in  Figure  24/  Font  2  is  a  Modified  fora  1  network  with  the 

e 

time  estimates  for  3  activities  lengthened  to  establish  one  dominant 
critical  path  through  the  network.  These  activities  are  also  noted  on 
the  network  diagram.  The  major  part  of  this  discussion  will  pertain  to 
the  network  referred  to  as  form  1. 

Although  the  approximation  formula  for  activity  mean  and  variance 
were  derived  from  the  Beta  distribution,  PERT  solutions  usually  consider 
activity  duration  to  be  normally  distributed,  as  well  as  project  duration. 
Simulation  of  this  network  (form  l)  was  accomplished  using  both  Normal  and 

a 

Beta  distributions  for  1000  trials. 

Histograms  of  the  resultant  distributions  for  the  two  simulation 
runs  are  given  in  Figures  25  and  26.  The  Chi  Square  Goodness  of  Fit 
Test  was  applied  comparing  the  simulated  resultant  distribution  with  the 
tabulated  Normal  Density  Function.  The  results  were  so  inconclusive, 
that  another  simulation  using  the  Beta  distribution  was  accomplished  for 
9999  trial  simulations.  A  histogram  of  the  resultant  distribution  of 
network  duration  for  the  9999  simulation  trials  is  given  in  Figui  '7. 

7 

Joseph  J.  Moder  and  Cecil  R.  Phillips,  Project  Management  With 
CPM  and  PERT  (New  York:  Reinhold  Publishing  Company) ,  p.  212. 
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TABIE  VI 

RESULTS  OF  THE  SIMULATION 
OF  THE  PERT  NETWORK 


Network 

Form  1 

Form  2 

Distribution 

Beta 

Beta 

Normal 

Beta 

Normal 

Number  of  Trials 

1000 

9999 

1000 

1000 

1000 

Mean 

49.6036 

49.6248 

49.5002 

58.6316 

58.9318 

Variance 

9.8428 

10.2629 

7.9085 

11.6842 

9.815^ 

X2  Value 

.92.3920 

1435.6596 

25.1858 

114.8033 

14.9355 

While  the  expected  values  (form  1  network)  are  all  essentially 

the  same,  the  variances  show  a  slight  increase  from  Normal  to  Beta 

2 

.distributions  for  1000  trials  to  the  Beta  9999  trials.  The  X  values 

for  both  simulation  runs  of  the  form  1  network  with  the  Beta  distribution 

are  highly  significant  indicating  that  normality  of  networks  containing 

2 

Beta  distributed  activities  could  not  be  assumed.  However,  the  X  value 
for  the  simulation  using  normally  distributed  activity  times  indicates 
acceptance  of  the  resultant  distributions  as  normal,  since 

X2  »  25.1858  <  X2_95j1?  =  27.587 

The  form  2  network  was  simulated  to  determine  if  a  network  with 


one  dominant  critical  path  would  improve  upon  the  results  of  the  pre¬ 
ceding  paragraph.  Three  activities  were  lengthened.  The  results  of  the 

1000  simulation  trials  utilizing  the  Beta  and  Normal  distributions  are 

2  "  , 

also  given  in  Table  VI.  While  the  X  value  for  the  Beta  distributed 
form  2  network  increased  to  a  value  more  highly  significant  than  that  of 


form 


it work , 


for  10C0 


trials 
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improved ,  as  indicated  by  a  decreased  in  the  value  of  X  . 

The  criticality  5ndex  for  all  simulation  runs  of  the  PERT  network 
(forms  1  and  2)  is  given  in  Table  VH,  indicating  the  relative  importance 
of  the  activities  within  the  network  contributing  to  network  realization. 
There  are  two  major  critical  paths  with  several  activities  in  common 
contributing  to  network  realization  for  the  form  1  network. 

Comparison,  of  the  information  given  by  the  criticality  index  and 
resultant  distributions  of  network  realizations,  reveals  that  both  forms 
of  the  network  containing  Beta  distributed  activity  times  produce  multi¬ 
modal  network  distributions.  The  resultant  distribution  of  the  form  1 
network,  where  there  is  more  than  one  critical  path  with  activities  in 
common,  more  closely  approaches  normality  than  the  resultant  distribution 
of  the  form  2  network,  where  only  one  critical  path  is  present. 

Another  factor  of  importance  concerns  use  of  the  PERT  assumptions 
in  determining  network  completion  time  defined  as  the  sum  of  expected 
activity  durations  on  the  critical  path  of  the  network.  The  critical 
path  is  the  path  of  maxxinum  sum  of  expected  activity  durations  in  the 
network  from  origin  to  terminal  nodes.  Based  on  this  assumption,  the 
expected  time  of  the  form  1  network  is  **7.65  and  variance  is  10.9. 
Comparison,  with  the  simulated  results,  indicates  the  PERT  estimate  is 
optimistic,  consistent  with  the  statement  made  by  Fulkerson  that  PERT 
produces  optimistic  estimates,  such  that  the  estimate  is  greater  than  the 
expected  value.® 

However,  the  PERT  assumption  of  normally  distributed  network  real- 

p 

D.  P. .  Fulkerson,  Exr? Cri  V ■- ?!  I-  ?T 

The  Rand  Corporation,  RM-3075-PE,  (Sania  Monica,  Karen,  p.  i. 


TABLE  VII 


CRITICALITY  INDEX  FOR  ALL 
SIMULATION  RUSS  OF  THE  PER?  NETWORK  ANALYSIS 

Criticality 


Activity 


For*  1  Network 


For*  2  Network 


I 

J 

Beta 

Beta* 

Hemal 

-  Beta 

Boreal 

1 

2 

0.00 

0.00 

0.00 

0.00 

0.00 

1 

3 

0.52 

0.51 

0.53 

0.81 

0.84 

1 

4 

0.48 

0.49 

0.46 

0.19 

0.16 

1 

5 

0.00 

0.00 

0.00 

0.00 

0.00 

2 

9 

0.00 

0.00 

0.00 

0.00 

0.00 

3 

8 

0.52 

0.51 

0.53 

0.81 

0.84 

4 

7 

0.48 

0.49 

0.46 

0.19 

0.16 

5 

6 

0.00 

0.00 

O.CO 

0.00 

0.00 

6 

11 

0.00 

0.00 

0.00 

0.00 

0.00 

7 

10 

0.48 

0.49 

0.46 

0.19 

0.16 

8 

10 

0.52 

0.51 

0.53 

0.81 

0.84 

9 

15 

0.00 

0.00 

0.0c 

0.00 

0.00 

10 

11 

O.65 

0.63 

C.?2 

1.00 

1.00 

10 

15 

0.35 

0.37 

0.28 

0.00 

0.00 

11 

12 

0.65 

0.63 

0.72 

1.00 

1.C0 

12 

13 

0.65 

0.63 

0.72 

1.00 

1.00 

13 

14 

0.65 

0.63 

0.72 

1.00 

1.00 

14 

25 

0.65 

0.63 

0.72 

1.00 

1.00 

15 

18 

0.26 

0.2? 

0.21 

0.00 

0.00 

15 

16 

o.ce 

0.10 

0.06 

0.00 

O.CO 

16 

17 

0.08 

0.10 

0.06 

0.00 

0.00 

17 

26 

0.08 

0.10 

0.06 

0.00 

0.00 

18 

29 

0.05 

0.05 

0.06 

0.00 

O.OC 

18 

21 

0.20 

0.22 

0.16 

0.00 

0.00 

18 

23 

0.00 

0.01 

0.00 

0.00 

0.00 

19 

20 

0.05 

0.05 

0.06 

0.00 

0.00 

20 

27 

0.05 

0.05 

0.06 

0.00 

0.00 

21 

22 

0.20 

0.22 

0.16 

0.00 

0.00 

22 

27 

0.20 

0.22 

0.16 

0.00 

0.00 

23 

24 

0.00 

0.01 

0.00 

0.00 

0.00 

24 

25 

0.00 

0.01 

0.00 

0.00 

0.00 

25 

2? 

0.66 

0.64 

0.72 

1.00 

1.00 

26 

27 

0.08 

0.10 

0.06 

0.00 

0,00 

27 

28 

1.00 

1.00 

1.00 

1.00 

1.00 

28 

29 

1.00 

1.00 

1.00 

1.00 

1.00 

♦This 

column  represents 

data  for 

the  9999  Simulation  trials 

of 

the  Beta  distributed  network.  All  other  data  is  for  1000  trials. 
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ization  times  utilizing  the  approximation  formula  for  mean  and  variance 
does  held  for  the  network  analyzed. 

This  chapter  has  presented  applications  of  the  GERT  Simulation 
model  through  simulation  of  various  GERT  networks.  The  PERT  assumptions 
were  investigated  through  analysis  of  GERT  networks  composed  of  all  AND 
nodes.  The  networks  utilized  to  test  probabilistic  GERT  applications 
required  approximately  one  minute  each  for  1000  simulation  trials  on  the 
CDC  3400  computer.  The  PERT  networks,  utilizing  the  Beta  probability 
distribution  required  approximately  4-  minutes  of  computer  time  for  1000 


simulation  trials. 


CHAPTER  V 


ANALYSIS  OF  AND  NODES 

The  results  of  a  GERT  network  analysis  are  the  probability  of  node 
realisation,  and  the  distribution  of  the  equivalent  time  parameter  for 
that  node  given  it  is  realized.  Extensive  research  has  been  conducted 
on  networks  containing  EXCLUSIVE-OR  logic  nodes,  while  conceptual  and 
computational  problems  exist  in  analysis  of  networks  containing  the  AND 
logic  node. 

Realization  of  the  AND  logic  node  occurs  when  all  the  activities 
inwardly-directed  to  the  AND  node  have  been  completed  or  realized.  Con- 
sideration  of  the  probability  of  node  realization  is  necessary  since 
GERT  networks  may  be  composed  of  probabilistic  and  deterministic  type 
logic  nodes.  Statements  about  the  associated  time  parameters  and  distri¬ 
butions  are  then  conditioned  on  the  probability  of  realizing  the  node. 
Should  the  network  consist  of  all  AND  nodes,  as  .in  PERT,  the  probability 
of  node  realization  is  one.  The  time  at  which  a  node  is  realized  is  the 
Maximum  or  the  durations  of  the  inwardly  directed  paths  to  that  node. 

The  purpose  of  this  chapter  is  to  investigate  general  applicabil- 
ity  of  the  merge  bias  correction  technique  proposed  by  Clark  in  analysis 
of  GERT  networks  containing  AND  nodes.  A  technique  is  required  to  com¬ 
pensate  for  erroneous  results  stemming  from  analysis  based  on  expected 
values,  as  in  PERT,  and  to  reduce  computational  difficulties  associated 

^Charles  E,  Clark,  "The  Greatest  of  a  Finite  Set  of  Random  Vari¬ 
ables,"  Oo,  Pes. ,  Vcl.  9>  No.  2,  I96I,  pp.  1^5-162. 
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with  the  maximization  operation.  Prt  ntation  of  this  chapter  is  by 
development  of  the  theory  and  through  application  on  two  sample  networks. 

The  merge  bias  correction  procedure  is  concerned  with  providing 
.  the  moments  of  the  maximum  of  random  variables  where  the  random  variables 
have  a  joint  probability  distribution.  Networks  composed  of  activities 
with  independent  distributed  variables  for  time  duration  provide  an  appli¬ 
cation  for  this  procedure. 

The  basic  assumption  of  this  procedure  is  that  the  independent 
activity  random  variables  w,  x,  y,  and  z  are  normally  distributed  with 
expected  value  u^  and  variance  G  and  G  Linear  correlation  may 
exist  between  any  two  paths,  such  as  x  and  y,  and  is  denoted  by  r(x,y)  or 
.  The  notation  0(x)  =  1/  If  2jf  exp(-x2/2),  and  i(x)  =  f  0(t)dt 
are  used  for  the  normal  distribution.  Also,  max(x,y)  is  denoted  by  x:y 

Let  v^  be  the  i^1  moment  about  zero  of  the  random  variable  max(x,y). 

The  following  equations  are  derived  by  Clark  and  are  utilized  in  solving 
for  the  maximum  of  random  variables. 


2 

a 


where  (T  ^  /  d*y, 


£  *  (ux  -  uy)/a 

V1  ®  *1  ft  )  +  u2  )  +  a  t  C  ft  ) 


4  i  . 


*  (u/  + 


')  )  +  (u22+  1  (-^  )  +  (1^  +  u2)  a  tip) 


The  linear  correlation  coefficient  is  given  approximately  as 
r(z  ,max(x,y))  =  s>x;j> 


! 
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Another  special  form  of  the  correlation  coefficient  is  also  use¬ 
ful. 

r(X  +  x,Y  +  y)  =  ^x+x,Y+y 

=  cJ/x  c/y  /°Xf y/^X+x  ^Y+y 

The  equations  given  above  my  be  utilized  in  obtaining  the  maxi¬ 
mum  of  any  number  of  random  variables  by  consideration  of  pairs  of  vari¬ 
ables  through  successive  application  of  the  following  procedure. 

Max  (  v,  x,  y,  z  )  =  Max  w,  max  (  x,  max  (y,z)  ) J 

At  this  point,  perhaps  a  simple  example  would  clarify  the  discus¬ 
sion.  Consider  the  network  shown  in  Figure  28  .  The  probability  of 

realizing  node  2  is  (.3)  (.*0  (.2)  (.8)  (.5)  =  0.0096  due  to  the  combin- 

2 

ation  of  AND  and  EXCLUSIVE- 0P.  nodes  as  constructed. 


FIGURE  28 

SAMPLE  NETWORK  UTILIZING  AMD  LOGIC  NODES 
The  time  of  realization  of  node  2,  T^,  given  it  is  realized  is  the 


a.  Alan  B.  Pritsker,  C-ERT :  Graphical 
Technique ,  the  FAIT)  Corpora t i : n ,  9" t 

Appor.dix  Ft. 


Evaluation  and  Review 

(S arra  Monica,  April, 


196o), 
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maximum  of  the  3  paths  leading  from  node  1  to  node  2  as  follows.  The 
notation  t^  represents  the  random  variable,  Time  of  duration  for  activity 

i. 

T,  =  Max  (t  f  t  ,  t  +  t  +  t  ,  t ,  +  t  5 
Z  a  d*  a  c  e  d  e 

?-  =  Max  (  t  +  t.  ,  t  +  Max  (  t  +  t  ,  t ,  )  ) 

2  a  o’  e  'a  c*d 

Traversing  the  network  for  T^,  assuming  node  1  has  a  mean  of  zero 
and  variance  of  one: 

=  t  ,  with  parairie  ‘-ers  Mean(A) ,  and 

Var(A). 


Node  A:  T. 


Node  C:  T, 


I 


-te(tA  +  T  ,  td) 


afc,d 

=0  ,  since  a+c  and  d  are 

parallel  paths. 

Vc 

=  u  +  u 

a  c 

wd 

*ttd 

2 

a 

fi  ■ 

=  d\2  +  cf  J 
d  a*e 

=  (  ud  -  Vc}  '  (ddZ  +  e/Jc  > 

V1  . 

=  ud  S  (,5  )  +  ua+o  it  (-{i)  +  a  t  ((3  ) 

v2 

-  <  “d2+  ^d2  >*  «j»>  +  <  Vo  +  d 

*.(^)+  Vo5  a  *  (r}- 

Mean  (C) 

=  vx 

Var  (C) 

2 

”  V2  "  V1 

a+-c 
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Node*  2:  T2  =  Max  (  TA  +  ^  ,  Tc  +  t0  ) 


Afb 


=  Mean(A)  + 


V 


d 


Afb 

/  2 
Cfe 


=  ua+  "b 
=  Hean(C)  +  u 

=  Var(A)  +  d 

=  Var(C)  +  cS' 


^Afb,C+e 

/“a.c 


"  f  A,'J  C'ktb  ^C+c 


a,d:afc 


C?d:afc  1 


The  correlation  coefficient.  O  .  =  0.  since  a  and  d  are  inde- 

r  a»d 

pendent  activities  in  the  network. 


y^a,afc 

^A,C 
P  Afb.DK 


a  v  a  /"a. a  ^  <1  ^afe 


=  <i  ^a  /l. 


afe 


=  o4fc  <4  5W^  !  ^  d: 


a+c 


Afb  C+e 


Thus  c  Iculation  of  a  new  a,  ,  v^,  and  v^,  for  node  2,  is 
possible  when  ^ Dfe  ^noWTi*  M®an(2)  and  Var(2)  follow  from  these 


calculations . 


) 
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The  equivalent  function  W  (s)  for  the  network  would  be. 

e 

w#(s)  =  p^2  exp(Mean(2)s  +  |  Var(2)s^) 
where  Is  the  probability  of  realizing  node  2,  calculated  at  the 
beginning  of  this  chapter. 

Another  procedure  would  be  to  add  pseudo  activities  in  the  network 

at  compensating  po'nts  to  visually  indicate  bias  corrections  affecting 

network  calculations  and  allow  conventional  (PERT)  network  solution 

3 

using  expected  times. J  The  pseudo  activities  would  in  a  sense  be  non¬ 
zero  dummy  activities  emanating  from  the  affected  node  with  mean  and 
variance  being  the  difference  between  merge  bias  and  conventional  calcu¬ 
lated  mean  and  variance. 

A  second,  more  complex,  network  with  numerical  values  is  offered 
to  further  illustrate  calculations  involved  and  simplifying  shortcuts 
which  may  be  used  to  decrease  the  computational  workload.  Since  the 
equations  utilized  are  the  same  as  the  first  sample  network,  only  a  brief 
coverage  will  be  given  the  actual  calculations.  The  GERT  Simulation  * 
Program  was  utilized  as  a  check  against  hand  calculation  and  offers  fur¬ 
ther  insight  into  solution  of  the  network.  Comparison  of  results  follows 
hand  calculations. 

The  network  is  as  shown  in  Figure  29;  where  the  activities  are 
a,  b,  c,  ***,k'f  node  numbers  are  1,  2,  3,  ***,6;  and  mean  and  variance 
are  denoted  as  (u,  <T  )  respectively.  The  probability  of  node  realization 
is  one,  so  all  activities  must  be  realized  fcr  network  completion. 


cp;i 


Joseph  J.  Moder  and  Cecil 
PERT  (liew  York:  Ro inhold 


R.  Phillips,  Project  bar a  cement  with 
Publishing  Co.,  19'-';),  n.  237. 
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Node  Is  Mean(l)  =  0  ,  Var(l)  =  0 

Node  2:  Kean(2)  =  4  ,  Var(2)  =  ls  since  activity 

independent  activity. 

Node  3^  =  Max  (  t  t  ,  t ,  ) 

3  a  c  d 


a  is  an 


Vo  =  7  .  »d  =  1  .ts'a  =  1 


0  =0  ,  since  a  +  c  and  d  are  independent  parallel 

!  «n"C  |d 

paths  with  no  common  element. 

2 


a 


=  2+1=3 


P 


=  -3.46 


=  6.997  ,  v2  =  50.9992 

Mean(3)  =  7 
Var(3)  =  2 

At  this  point  a  simplifying  aspect  is  noted,  in  that  Mean(3)  and 
Var(3)  are  at  most  l/lO  per  cent  off  the  conventional  calculations  using 
expected  values.  A  rule  given  by  Moder  and  Phillips  applies  in  this  case:1 

If  the  difference  between  the  expected  complete  times  of  the  two 
merging  activities  being  considered  is  greater  than  twice  the  larger 
of  their  respective  standard  deviations,  then  the  bias  correction 
will  be  small,  less  than  a  few  per  cent,  and  can  be  ignored. 


Node  4:  T., 


u 


2fb 


=  Max  (T^  +  tQ ,  T2  +  t^) 


=  6  *  u3*-e  =  10  »  ^  2+b  =  2  »  *  >e 


&  2+b‘,>e  =  *41 


a 


=  3  ,  /?  =  2.31 


Ibid. ,  p.  239. 
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II 

% 

9.99  , 

v2  =  102.92 

Moan  (4)  = 

9.99,  Var(4)  =  2.92 

Here  the  rule  also  holds,  but 

to  a  lesser  effect,  with  the  bias 

correction  at  .01  and  .08  for  the  mean  and  variance  respectively. 

Bode  5: 

Max  (  T4 

+  tg,  t2  +  tf  ) 

U2ff 

9,  Vg  : 

■  U.  <*\t  - 2-  tS'Vg^3-92 

f  z* {Me,  = 

.439 

2 

a  = 

2.93,  { 

?  =  1.17 

it 

11.1  , 

v2  =  126.8 

Mean(5)  = 

11.1,  V 

t  (5)  =  4.78 

Node  6:  = 

Max  (  T4 

+  V  T5  +  tk  ) 

Vh 

13.99, 

u5fk  =  *  ^4+h  " 

^  5tk 

=  5.78 

r 4hb,5hk  = 

.286 

2 

a  = 

7.0 

?  =  .336 

vi 

14.68, 

v2  s  218.79 

Mean(6)  = 

14.68, 

Var(6)  =  3.29 

• 

The  sanple  network  was  also  simulated  using  the  GERT  Simulation 

Program.  All  activities  we 

re  normally  distributed  with  mean  and  variance 

as  indicated  in  Figure  ?9  on  cage  9-  . 


The  number  of  s inula t ion.  trials 
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was  1000.  A  favorable  comparison  of  the  hand  calculations  and  simulation 


results  is  given  as  follows 

for  nodes 

4  and  6. 

Node 

4 

6 

Simulation 

Mean 

10.048? 

14.3505 

Variance 

3.0128 

3.6519 

Hand  calculation 

Mean 

9.99 

14.68 

Variance 

2.92 

3.29 

The  criticality  index  indicates  one  path  is  predominant,  approx¬ 
imately  64  per  cent  of  the  trials,  through  the  network  from  nodes  1,  2, 

3,  4,  to  node  6.  A  minor  critical  path  from  nodes  1,  2,  3,  4,  5>  to  node 
6,  occurs  approximately  31  per  cent  of  the  time.  A  histogram  of  node 
realization  times  for  nodes  4  and  6  is  given  in  Figure  30. 

The  assumption  of  normality  for  non-normal  activities  introduces 

error  into  the  calculation  and  is  discussed  at  length  in  the  referenced 

5 

article  by  Clark.  Approximation  of  the  first  two  moments  and  v^  of 
Max(x,y),  a  non-normal  distribution,  with  the  first  two  moments  of  a  nor¬ 
mal  distribution  does  not  produce  too  great  an  error  when  the  expected 
values  and  variances  the  same,  for  both  distributions. 

An  example  is  given  where  the  greatest  of  500  standard  normal  var¬ 
iables,  which  is  not  normally  distributed,  is  approximated  by  a  normal 
variable  with  negligible  error.  Approximations  of  exponentially  and 
uni formally  distributed  random  variables  by  normally  distributed  vari¬ 
ables  are  also  exhibited  resulting  in  a  larger,  but  tolerable  error. 

5Clark,  on.  cit.,  n.  151. 
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The  procedure  of  merge  bias  correction  does  provide  a  means  for 
overcoming  errors  inherent  in  conventional  network  calculations  based 
on  expected  values  as  in  PERT.  The  calculations  are  tedious  for  hand 
calculation,  but  adaptable  to  computer  solution. 

The  procedure  is  applicable  to  GERT  network  analysis,  where  the 
sub-re tworks  containing  AND  nodes  may  be  reduced  to  equivalent  functions. 
Should  the  complexity  of  the  network  between  AND  nodes  prevent  network 
reduction,  addition  of  the  pseudo  activities  to  the  network  would  facil¬ 
itate  solution  by  the  maximization  operation  using  expected  values. 
However,  general  applicability  of  merge  bias  correction  cannot  be  claimed 
for  GERT  networks,  as  this  procedure  is  mainly  problem  oriented. 


CHAPTER  VI 


SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 
I.  SUMMARY 

The  purpose  of  this  research  was  to:  (1)  develop  a  general  purpose 
GERT  Simulator  allowing  analysis  of  GERT  networks  through  the  technique 
of  Simulation,  (2)  investigate  resultant  distributions  of  GERT  networks 
containing  all  AND  logic  nodes  to  determine  validity  of  the  PERT  normal¬ 
ity  assumptions ,  and  (3)  investigate  analysis  of  AND  logic  nodes  through 
analytical  and  simulation  methods  for  distribution  of  the  equivalent  time 
parameter . 

The  goal  to  develop  a  generalized  technique  for  analysis  of  complex 
systems  portrayed  by  stochastic  networks  has  not  been  fully  realized.  The 
general  purpose  GERT  Simulation  Program  was  needed  to  assist  the  research 
endeavor  and  reduce  the  computational  complexity  in  solution  of  large 
networks . 

The  Simulation  Model,  was  designed  to  facilitate  solution  of  GERT 
networks,  allowing  as  much  freedom  as  possible  in  depicting  networks  for 
analysis.  The  program  is  user-oriented,  as  input  data  is  held  to  a  min¬ 
imum,  requiring  only  one  input  card  for  each  activity  in  the  network. 

Five  probability  density  functions  are  available  for  use  in  describing 
the  distribution  of  activity  durations.  Program  output  includes  a  criti¬ 
cality  index  on  each  activity  as  the  relative  frequency  an  activity 
appears  on  the  critical  path.  The  mean  and  variance  of  network  completion 
is  also  provided  as  output,  along  with  analysis  of  specified  nodes,  con¬ 
sisting  of  mean,  variance,  probability  of  realization  and  histograms  of 
node  realization  times.  A  Chi-Square  Goodness  of  Fii  Te-1'  mv  bo  a*  olio i 


to  tbe  distribution  of  node  realization  times  and  printed  as  output. 
Various  GERT  networks  were  simulated  to  provide  examples  of  the  use  of 
the  Simulation  Model  and  verify  results. 

An  investigation  into  the  PERT  normality  assumptions  was  conducted 
through  use  of  the  Simulation  Model.  The  Beta  and  Normal  probability 
distributions  were  utilized  in  simulating  a  sample  network.  Resultant 
distributions  of  network  realization  times  were  analyzed  to  determine  if 
normality  could  be  assumed. 

Analysis  of  AND  logic  nodes  was  conducted  using  simulation  and 
analytical  techniques.  The  merge  bias  correction  procedure  was  applied 
to  sample  networks  to  determine  if  the  procedure  was  applicable  to  solu¬ 
tion  of  GERT  networks  containing  AND  nodes. 

II.  CONCLUSION 

The  GERT  Simulation  Model  did  satisfy  all  requirements  set  forth 
in  Chapter  I.  Simulation  of  sample  GERT  networks  provided  logical  and 
conclusive  results  as  established  by  appropriate  statistical  tests. 

The  normality  assumptions  utilizing  approximation  formulas  for 
mean  and  variance  of  activity  duration  in  solution  of  PERT  networks  was 
verified  by  analysis  of  resultant  network  distributions.  Strict  appli¬ 
cation  of  the  Beta  probability  distribution  for  activity  duration  times 
did  not  produce  normally  distributed  network  realization  times  for  either 
of  the  network  configurations  considered. 

•The  merge  bias  correction  procedure  utilized  to  produce  better 
approximations  of  network  realization  times  than  conventional  methods 
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although  providing  accurate  estimates  of  distribution  parameters. 

However,  the  normality  assumption  of  activity  distribution  ccuid  intro¬ 
duce  significant  errors  in  estimation.  This  procedure  is  therefore 
problem-oriented  and  not  applicable  to  general  solution  of  GERT  networks. 

HI.  RECOMMENDATIONS 

The  use  of  the  Simulation  Model  will  facilitate  further  research 
into  GERT  networks.  However,  further  applications,  other  than  those 
presented  in  this  paper,  are  necessary  to  assure  satisfactory  results 
in  tho  general  case,  and  are  therefore  recommended. 

The  application  of  the  merge  bias  technique  to  general  GERT  net¬ 
work  solution  for  networks  containing  AND  logic  nooes  is  computationally 
difficult.  Further  investigation  into  this  procedure  and  the  underlying 
normality  assumptions  is  recommended.  A  computer  program  to  solve  PERT 
networks  wioh  this  procedure  would  provide  resultant  statistics  for  com¬ 
parison  with  simulated  and  analytical  results. 
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